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ABSTRACT 

The statistical technique of analyzing data 
collected at reqular intervals to tevcal periodic components of the 
data is described by reference to actual records. The data chosen for 
illustration include tide height in a river; biochemical oxvaen 
demand and dissolved oxygen in the same river; discharged salt into a 
river system and its relation to the hardness of the water; and the 
discharge of sewage into a river. Annual, biannual, and other cvcles 
are demonstrated in the data. The limitations of the technique, the 
sources of error resulting from methods of data collection ant choice 
of parameters for use in the statistical analysis are discussed. It 
is emphasized that interpretation of the results must be based on an 
understanding of the natural process, "’'he book is primarily intended 
for use in mathematical modeling of environmental parameters in 
studies of water pollution of rivers and estuaries. There is a short, 
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PREFACE 

T his volume is designed to offend almost everyone engaged in the control 
of water pollution. Sanitary engineers will be offended because there is 
too much mathematics; mathematicians, statisticians, and systems analysts 
will object because the mathematics is over-simplified; administrators and 
managers will point out that there is no panacea here for management 
decisions; oceanographers will be shocked because the author has felt it 
necessary to explain a technique which has been common practice for 
many years in oceanography. 

All of these criticisms are just. Spectral analysis has been used for many 
years in oceanography, meteorology, electrical engineering, medicine, and 
several other fields. It is a method of data analysis and, as such, is a con- 
tributor toward the technical foundation upon which management can 
make a decision, but it provides no magic mathematical model to assist in 
the process. The mathematics and statistics of what is essentially a mathe- 
matical and statistical tool arc indeed over-simplified; but the engineering 
rationale is also over-simplified t > make it comprehensible to mathema- 
ticians and systems analysts. 

In short, the author has attempted the dangerous task of bridging fields 
and trying to explain one field to another. The concepts and techniques of 
computation presented here are merely the author's interpretation of what 
others have developed — for these ideas he takes no credit. The interpre- 
tations of spectral results in terms of sanitary engineering phenomena are 
entirety the author’s ideas — for these he takes all the blame. 

The past decade has seen great changes in the basic approach to the 
study of water pollution and water resources problems. The advent of 
high-speed digital computers has resulted in a spate of mathematical 
models based on the use of the rather sparse data obtained by manual 
techniques. At the same time advances in instrumentation, recording, and 
telemetry have made remote sensing stations operational in environmental 
monitoring, but the use of data from them has been minimal in water 
pollution control. 

It is the gap between environmental monitoring and mathematical mod- 
eling that the author has attempted to bridge by pointing out the applica- 
tion of spectral analysis in the interpretation of time series data, both in 
estuaries and in free-flowing streams. 

Spectral analysis as a technique emerged from the work of La Place 
and Fourier near the beginning of the nineteenth century, and much of 
the mathematical theory underlying its practical application was developed 
in the work of 0. I. Taylor, Norbert Wiener, and S. O. Rice. With the 
advent of high-speed digital computers and the development of automatic 
sampling techniques its practical use in many fields became feasible. The 
author offers his own experiences in the use of spectral analysis in water 
pollution problems as a first attempt in the use of this technique in the 
study of water pollution in rivers and estuaries. 

T. A. WASTLER 
WxsmNoroN, D. C 
DtctMitA IMS 
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THE CONCEPT OF SPECTRAL ANALYSIS 



O ne of the classic experiments in the study of optics is to pass a beam 
of sunlight through a triangular prism and observe the rainbow of 
colors into which the beam of white light is split. This experiment serves as 
a useful physical analogy to the operation of spectral analysis on a body 
of data. 

In Figure la is presented a schematic diagram of the results of passing 
a beam of sunlight through a triangular prism. The beam of sunlight is 
split into a spectrum of colors ordered according to their respective wave 
lengths (or frequencies). If this spectrum of colors is allowed to strike a 
battery of light-sensitive cells, the intensity of the light of each color can 
be measured; these results can then be plotted as a “light intensity spec- 
trum, which might look like that in Figure la if each of the six principal 
colors had the same intensity. 

If a similar experiment is performed with a light beam that Is composed 
of only three of these colors (present at different intensities), the resulting 
light frequency spectrum and light intensity spectrum might resemble 
those in Figure lb. 

This experiment demonstrates the resolution of a complex physical 
phenomenon into a group of simpler phenomena that may be easier to 
examine from both theoretical and practical viewpoints. 

The effect of using spectral analysis on a record of observed field data 
is directly analogous to the effect of the prism on the light beam. This 
analogous effect is presented schematically In Figure Ic. The actual tech- 
nique of spectral analysis is discussed later; at this time the significance 
of the result of the computation is of concern. 

In the prism experiment the relative intensities of the resolved light fre- 
quencies can be observed and studied. Spectral analysis of a record of 
observations results in a sorting of the total “variance" of the record into 
its component frequencies. The variance of a data record is therefore 
analogous to lire intensity of the light beam. 

The vorianct is defined as the sum of the squares of the deviations from 
the mean divided by one less ’tvn the number of observations. This is the 
definition of variance as it it no mally used as a descriptive statistic of a 
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Figure 1. Physical Analogy to Spactrat AnatyiH. J 
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body of data. Conceptually, the variance is a measure of the dispersion 
of reservations about the mean value. In the statistical interpretation of 
data, it is ordinarily regarded that this dispersion of values about the mean 
is due to chance. 

The statistical manipulation of time-series data by spectral analysis 
results in the computation of those parts of the variance of a record that 
recur at constant time intervals as well as the part that it random (non- 
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recurring) in character. As the light beam may be resolved by the prism 
into its component colors of different intensity, so is the variance of a 
time-series record resolved into its component parts by spectral analysis. 

The estimates of variance for each frequency resolved in the spectral 
analysis form the “power spectrum” of the record front which the compu- 
tations are made. (The term “variance spectrum" would be more accurate; 
but in the pioneering work in spectral techniques done in communications 
engineering the term "power," which is closely related to record variance 
in that frame of reference, came into common usage.) The power spec- 
trum computation is the fundamental operation of data reduction and 
interpretation by means of spectral analysis; the many other computations 
that can be made in spectral analysis are all based firmly upon the con- 
cept and calculation of individual power spectra. It may be stated that the 
computation of individual power spectra bears about the same relation- 
ship to spectral analysis as differentiation and integration bear to the Cal- 
culus. 

The interpretation of variance as a statistic descriptive of both the ran- 
dom and nonrandom characteristics of a time-series data record is most 
important in understanding the significance of spectral results. In the usual 
type of statistical analysis, variance is conceptually regarded as a measure 
of the random dispersion of the observations from their mean value. In 
many cases this is true; but that this is not a necessary condition for the 
existence of a variance can be demonstrated with the aid of Figure 2, in 
which segments of three hypothetical records and the corresponding power 
spectra are presented. 

Figure 2a shows a record that has a constant value, i.e., all values are 
equal to the mean. Since there are no deviations from the mean, the vari- 
ance is rero and the power spectrum is mo at all frequencies. 

Figure 2b shows a record that forms a sloping straight line. The seg- 
ment of the record shown in this figure has a mean of 3.55 and a variance 
of 0.69. It is apparent that none of this variance results from a “random" 
dispersion about the mean, but that the variance is the result of a secular 
(time-dependent) trend in the record. If spectral analysis were done on 
the record of which this segment is a pari, all of the variance (or "power") 
would be concentrated in the zero-frequency spectral estimate as shown 
in Figure 2b. The zero-frequency spectral estimate includes all of the 
record variance that does not recur during the length of the record used 
in the analysis. It therefore includes (I) any truly random fluctuations in 
the record, (2) any linear trends in the record, and (3) any periodic 
components in the record that are of so low a frequency that they appear 
as linear trends in the record. For example, spectral analysis of the seg- 
ment “A" in Figure 2c would result in a power spectrum similar to that 
in Figure 2b, simply because the record length is not great enough to 
resolve the periodic fluctuation exhibited in Figure 2c. 

The record of which Figure 2c is a small part also exhibits a variance. 
The mean of this record is 4.0, and the variance is entirely the result of 
the sinusoidal fluctuation about the mean. Spectral analysis of this record, 
which should be about 10 times longer than the segment presented in 
Figure 2c, would result in a power spectrum In which all of the variance 
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is concentrated at the frequency designated f* which corresponds to the 
wave length of the sine wave in the data record. 

If the data record were a combination of Figures 2b fc/id 2c, the power 
spectrum would be a combination of the spectra in these figures, Le„ there 
would be components at the zero frequency and at f 2 . It is this charac- 
teristic that makes spectral analysis such a useful tool in analyzing records 
that represent complex phenomena in natural systems. For example, the 
diurnal effects of photosynthetic activity could be separated from the 
longer-period effects of waste loads and river discharges by spectral anal- 
ysis of the stream dissolved oxygen (DO) record. 

The power spectra in Figures 1 and 2 are presented in bar graphs to 
emphasize two characteristics of spectral results: (1) estimates of the 
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variance at several discrete frequencies arc obtained from the analysis and 
(2) each variance estimate represents the variance concentrated in a band 
around the nominal frequency of each variance. This bar graph represen- 
tation is not the usual way in which spectra are presented; in the remain- 
ing figures discussed in this paper, the conventional point-and-line repre- 
sentation is used. 

Like any other method of data analysis, this method has its limitations 
and disadvantages. Three major requirements for the record may be 
regarded as limiting. First, the record must be fairly long — generally hav- 
ing over 100 sequential measurements. Second, there must be no missing 
data — if measurements are missing, suitable values must be interpolated 
before spectral analysis is attempted. Third, the mathematical procedures 
require so much computation that the use of a high-speed digital com- 
puter is essential for most analyses. These restrictions are discussed in 
more detail. 
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THE TECHNIQUES OF SPECTRAL ANALYSIS 



T he approach used here is to present the technique of spectral analysis 
by following through the steps in the actual spectral analysis of a record. 
For the benefit of those who wish to explore the mathematical basis for 
spectral analysis computations, pertinent references are presented in the 
bibliography. 

The record chosen for detailed examination here is a record of water 
level obtained from a U.S. Geological Survey station in the Potomac 
Estuary near Washington, D. C. This record was chosen because it exhi- 
bits a simple periodicity with very little random interference (“noise”), 
and because the data were obtained as a continuous recording so that a 
wide choice of sampling intervals was possible. 

A portion of this record is presented in Figure 3. A visual examination 
of this record shows that there is a dominant period of about 12 hours 
and that there is some long-period change. 

Although the entire computation can be carried out on high-speed digi- 
tal computers, with available programs, the individual steps are presented 
here to illustrate the technique. Only the initial steps in the data prepara- 
tion need be carried out by hand or on semi-automatic equipment. 

Step 1 . A sampling interval of four hours was chosen for this particular 
analysis, and the record was read at this interval for a total of 145 read- 
ings. (Considerations governing the number of points read and the sam- 
pling interval are discussed later in detail.) The starting point was arbi- 
trary. The numbers obtained, in order, are 
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Figure 3. Portion of Tide Height Record in 
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Step 2 , The autocorrelation function of these numbers is then formed. 
This is a very large name for a very simple, but very useful, process. Each 
number in the record is multiplied by another number in the record, and 
from the mean of the sum is subtracted the square of the arithmetic mean 
of the entire record. If the individual products are formed by multiplying 
each number by itself, the result is called the “autocorrelation at lae 0,” 
which is purely and imply the “variance” as ordinarily defined in statis- 
tics. If the individual products are formed by multiplying each number by 
the number that follows it in the sequence, the result is called the “auto- 
correlation at lag 1.” The autocorrelations computed to 1? lags for the 
record being analyzed here are 
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C 3 a Autocorrelation at lag 2 
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C 8 s Autocorrelation at lag 3 
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The remaining autocorrelations are computed similarly and have the 
values 

C 4 = 0.154 C f = 0.0476 C 10 = 0.252 

C 8 = 0.735 C s = -0.917 C n = -0.998 

C e = 1.10 Co = 0.921 C 1S = 0.798 

This operation may be expressed mathematically as 

1 n ~ r T i ” l! 

Cr = — E ~ [ n f ** ] , (1) 

where 

C, = autocorrelation at lag r, 
x t = record value at /, 

t = 0, 1, 2 ... n = sequential index of values, 

r = 0, 1, 2 ... m = lag number, 

m = total number of lags. 



Step 3. The Fourier cosine transform for each autocorrelation is next 
computed. This serves to smooth out some of the wide fluctuations in the 
autocorrelations and consists of applying a cosinusoidal weighting factor 
to each autocorrelation calculated in the preceding step. This operation 
can be expressed mathematically as 
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V, 



( 2 ) 



where 




[ m— 1 

Co4~^»n cos r«— 1-2 2 » 

9 =i m j 



F r = Fourier cosine transform of the autocorrelation at lag r, 
q = lag number, having values between 1 and m— 1, 

A = a constant, 

k “ 1 for r = 1, 2 . , , m— 1, 

A; = V4 for r = 0, r = m, 



and the other letters have the definitions previously given. 

The Fourier cosine transforms calculated for each autocorrelation of the 
tide height record are 

Vo = (2)712) [ +0 . 7 98 cos(0)* ] 

q 1 

= ^ £ 1.32+2(1) (—.312— .553+1.20— 0.154— 0.735 

+ 1.10+0.0476-0.917+0.921+0.252-0.998) 

+0.798(1) J 
= 0.0765 

Vi = J2 [ ^ .32+ (0.798 )cos(l )jt ] + -| [-(0.312)cos f 1 
-(0.553)cos^^+(1.20)cos ^^-(0.154) 
eos^^y^ — (0.735) COS +( 1.10) cos 

+ (0.0476 ) cos (7) 1 ( 2 1)ir - (0.9 1 7 ) cos - + (0.921 ) 

col J^ +(a252)co( n^_ (a9 98)co.t^ ] 

V, = 0.101 



Similarly, the Fourier cosine transforms for the remaining autocorrela- 
tions can be computed: 



V t = -0.0332 
Vo = 0.0570 

Vi = -0.0443 



Vi = 0.995 

Vi = -0.224 
K 10 = 0.146 
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V 6 = 0.0943 

Ve ~ -0.136 
Vi = 0.353 


V n “ -0.118 
V it = 0.0561 


Step 4. The final step in the spectral analysis of a single record is an- 
other weighting operation that counteracts some distortion of the spec- 
trum resulting from the small sample size. 

This step can be expressed mathematically as 


U 0 = 0.54 [V 0 + Vil 


(3) 


U r = 0.2'iV r -i-\-0.54V r -\-0.23V r +i, 


(4) 


for r — 1,2,3, . . . m— 1 
U„ = 0.54 *' m _ 1 +0.54 V m 


(5) 



where t/ 0 , U n Um are the power spectrum estimates corresponding to the 





Flflurf 4. Spectrum of a TkJal HeJofit Record 



12 



respective lags, and the remaining symbols have the meanings previously 
assigned. 

The power spectrum estimates for the tide gage record are 



u o 


= (0.54) (0.0765) -+-(0.54) (0.101 ) = 
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= 0.0183 






U 3 


= 0.0130 


Us = 


0.567 


V 4 


= 0.0109 


U 9 = 


0.142 


u s 


= 0.00947 


Uio = 


0.000446 


c/ 8 


= 0.0295 


U n = 


-0.0173 


c/t 


= 0.388 


C/,2 = 


-0.0334 



Each of these spectral estimates represents the part of the total record 
variance that is estimated to occur with a certain periodicity. The period 
corresponding to each lag is determined from the lag number and the 
sampling interval by this relation: 

Tr = ~~ , ‘ (6) 



where 

T r = period corresponding to lag r, 

Ar = sampling interval. 

For the record analyzed here, m — 12 and Ar = 4 hours. The periods 
corresponding to the lag numbers are these: 



Number 


Period (hours) 


0 


00 


1 


96 


2 


48 


3 


32 


4 


24 


5 


19.2 


6 


16 


7 


13.7 


8 


12 


9 


10.7 


10 


9.6 


11 


8.7 


12 


8 



The spectral estimates are plotted as functions of period in Figure 4. 

From these results is can be seen that about 79 percent of the observed 
total variance in water level can be attributed to periodicities of around 
12 hours, that there are very small diurnal effects in the tide at this point, 
and that long-period changes account for about 12 percent of the ob- 
served variance. 

The effects of using different sampling intervals and computing to dif- 
ferent numbers of lags may be examined by considering the results of 
further spectral computation with this record, Figures 4, 5, and 6 illus- 
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(rate spectra obtained from the same record, with variations only in the 
sampling interval and in the number of lags used in calculation. The total 
record length in each case was 576 hours. The spectral estimates from 
which Figure 4 was plotted were obtained by computation to 12 lags from 
values read at 4-hour intervals (145 points). Figure 5 was obtained from 
computation to 12 lags from values read at hourly intervals (577 points). 
Figure 6 was obtained from computation to 24 lags from values read at 
hourly intervals. The abscissal scales in Figures 5b and 6b are 50 times 
greater than those in Figures 5a and 6a; this magnification shows short- 
period effects more clearly. 

While each of these figures shows the dominance of approximately 
semidiurnal periodicities in this record, the effects of the different sampling 
intervals and the different numbers of lags can be seen also. 

Comparison of Figures 4 and 5 shows the effects of a change in sam- 
pling interval on the spectrum of the tide height record. First, the use of a 
smaller sampling interval for a given record length increases the number 
of measurements used in the analysis and therefore increases the number 
of degrees of freedom upon which each estimate is based. This results in 
the smaller confidence band shown in Figure 5. Second, the change of 
sampling interval from 4 hours to 1 hour permits the resolution of com- 
ponents with periods as short as 2 hours in the latter case instead of the 
8 hours possible in the former. However, the use of the same number of 
lags with the shorter sampling interval does not permit as high a degree 
of resolution of long-period phenomena as was obtained with the longer 
sampling interval. In fact, the spectrum calculation leading to Figure 5 
does not permit an estimate of the diurnal and longer-period effects be- 
cause these are effectively masked by the dominant 12-hour component. 

Comparison of Figure 5 with Figure 6 shows the effect of increasing 
the number of lags without changing the sampling interval. The resolving 
power is increased, i.e., estimates are made for a greater number of spec- 
tral bands, but with a slight loss of confidence in the estimates, as is indi- 
cated by the wider confidence band in Figure 6 than in Figure 5. In this 
particular case, however, the increased resolving power demonstrates the 
existence of a small but statistically significant overtide (a harmonic of the 
semidiurnal tidal component) with a period of about 6 hours. This over- 
tide is not shown in Figure 5b because the dominating semidiurnal band 
has overlapped the 6-hour-period band sufficiently to mask the very small 
overtide. In Figure 6b the resolving power of the 24 lags used in compu- 
tation produces band-widths sufficiently narrow to prevent ovr rlap of the 
12- and 6-hour periods. The existence of higher harmonics is also shown 
in both Figures 5b and 6b. The former figure shows significant fourth and 
fifth harmonic overtides, whereas the latter shows these and several other 
short-period components. 

The results of the spectral analysis of this record may be used to point 
out several characteristics of this technique. 

First, each of the spectral values obtained represents an estimate of the 
variance over a range of periods in the vicinity of the nominal value. The 
range of periods for which each estimate is computed is determined by the 
sampling interval and the number of lags used in computation. This range 
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is called the "equivalent width" (IV,) of the spectral band and can be cal- 
culated in terms of frequency by 



The spectral estimate reported is an average value for all periods in the 
band over the range IV r . For example, for the results presented in Figure 
4, the spectral estimate reported as corresponding to a 12-hour period is 
actually an average for the range of periods from 10.7 hours to 13.7 
hours. The overlapping of spectral estimates is illustrated in Figure 5, and 
Figure 6 demonstrates how this overlapping can be reduced by increasing 
the number of lags, thereby reducing the equivalent width of the band for 
each spectral estimate. 

Second, the precision of each estimate is a function of the total number 
of samples and the number of lags used in computation. A method for 
estimating the number of degrees of freedom for each estimate and for 
establishing the confidence intervals has been presented by Bl? ckman and 
Tukey. In this method, the process being measured is regarded as Gauss- 
ian and the degrees of freedom and confidence intervals are based on a 
Chi-square distribution. The 80 percent confidence bands indicated on 
Figures 4, 5, and 6 have been determined by this method. It has been pre- 
viously noted that a decrease in the total number of samples or an increase 
in the number of lags causes a widening of the confidence band, which 
represents a loss of precision in each spectral estimate. 

Third, if the sampling interval is not small enough to permit resolution 
of the shortest periods that contribute significant variance to the record, 
the short-period variance in the record is not lost but Is reported at har- 
monics of the true period. For example, if this record, as analyzed in 
Figure 4, contained significant variance with a period of 4 hours, this may 
be shown in the computations as part of the spectral estimate for the 8- 
hour period. This occurrence is known as "aliasing" or "folding." In 
Figure 4, the appearance of a significantly large estimate (in terms of 
absolute magnitude) at the shortest period computed in the analysis indi- 
cates that there is probably some aliasing of the short-period variance. 
The use of a smaller sampling interval, as shown in Figure 5, eliminated 
the aliasing in this spectrum. 

Fourth, the appearance of an occasional negative spectral estimate is 
an artifact of the computation and results from the use of a finite record 
length for spectral computations. If the record were infinitely long, there 
would be no negative spcc'ral estimates; but, when an estimate is close 
to zero in magnitude, it may appear with either a positive or negative sign. 
This result is interpreted only as a very small quantity of variance, and 
no physical significance is attributed to the negative sign. 
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INTERPRETATION OF SPECTRA 

T o illustrate how the computed spectra of a set of observations can be 
used to gain an insight into the structure of a river or estuarine system, 
some data obtained in a field survey of the Potomac Estuary near Wash- 
ington, D. C., are presented here as a basis for discussion. Since the pur- 
pose of this discussion is to present a conceptual picture of how spectral 
results may be interpreted, the quantitative results obtained from this 
particular study are not presented at this time. 

In Figures 7 and 8 are presented time series of dissolved oxygen (DO) 
and 5-day biochemical oxygen demand (BOD) measurements obtained 
at six stations in this estuary. In Figures 9 and 10 are the spectra com- 
puted from these records. Each record contained 145 points obtained at 
4-hour intervals; computation was carried out to 12 lags. Figure 11 is a 
schematic representation of the station locations. 

Each spectrum shown in Figures 9 and 10 consists of three major 
effects: long-period, diurnal (24-hour), and semidiurnal (12-hour). The 
estimates adjacent to the long-period, 24-hour, and 12-hour periods are 
affected by the strength of the dominating period and the overlapping of 
spectral bands, as was discussed earlier in relation to the tide gage record. 
The significance of each of these effects may be considered separately. 

The effects reported in the computation as '‘long-period’ 1 or of "infinite” 
period include all components whose periods are too long to be resolved 
in the computation. Since in these computations there is overlapping of 
the 96-hour and "infinite" periods, this means that all components with 
periods longer than 48 hours are regarded as long-period components. 
The use of the infinity symbol in the figures is merely a convenience in 
plotting and should be interpreted as meaning "long-period." 

The components reported as long-period may have one or more of 
these physical interpretations: 

I . The record may be affected by the existence of periodic components 
that art too long to be resolved by the length of record available. 
For example, a record of air temperature for a length of several 
months might show a long-period effect that, from a record of 
several years duration, would appear as an annual cycle. 
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Figure 7. OiaioNvd Oxygsn Records Obtained 'm th# Potomsc Estuary, August, 1959. 
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Figurv & Biochemical Oxygan Demand Record Ofttt fa nd In the Po tt nu c Efoivy, August. 1969, 
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Figure 11. ScMmttte RtpreMftttttoft of Potomac Ertuary Sampling Station Locations. 

2. There may be a secular trend in the record that l< fundamentally 
aperiodic in nature. For example, a record of world population 
over several centuries might show such a trend. 

3. Random sampling, reading, and analytical errors appear as long- 
period effects. A constant bias in the data affects the mean of the 
record only, and since spectra are computed from deviations from 
the mean, such a bias would not affect the spectral results. 

From th » computed spectra alone, it is not possible to differentiate 
among random errors, aperiodic effects, and long-period effects. In a 
practical tense it is often quite adequate to regard the long-period spectral 
estimates as secular trends for the available record length and to regard 
any random error as constant for all records of the same measurement. 
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The DO spectra presented in Figure 9 exhibit long-period effects that 
generally decrease from the upstream to the downstream stations. If a 
constant random error is assumed, the existence of an effect that is a func- 
tion of distance from the head of tidewater may be postulated, and it 
might be deduced that this effect is related to the river discharge entering 
the estuary above Station 3. Since the parameter being measured is DO, 
this long-period effect might be interpreted as a measure of the amount 
of DO advected to the system in the river discharge or as an effect of river 
discharge on the reaeration characteristics of the system, or both. Con- 
siderable caution is in order if these results are to be interpreted in this 
manner; the change in the long-period estimate from Station 3 to Station 
8 is barely significant at the 80 percent level, and the differences in this 
estimate between successive stations are not significant at the 80 percent 
level. The existence of a definite trend does indicate, however, that the 
interpretation given is reasonable and affords a basis for the examination 
cf other results in terms of this hypothesis. 

The BOD spectra in Figure 10 present an aspect somewhat different 
from the DO spectra at the same stations. The long-period effects at four 
of the six stations have very similar values. 

Station 7 exhibits a long-period component that is somewhat higher 
than these but still lies within the 80 percent confidence band of the four 
stations. Station 3 exhibits a long-period effect that is significantly less 
than the corresponding effects at the other five stations. These components 
at the five downstream stations can be interpreted as the result of a secu- 
lar trend in the major waste discharge to the system as well as random 
error in the determination of BOD. 

Comparison of the long-period effects at Stations 3 and 4 with those of 
stations in the immediate vicinity of the waste outfall suggests that the 
limits of a tidal excursion upstream from the outfall may actually lie be- 
tween Stations 3 and 4 instead of downstream from Station 4, as the 
calculations based on the tidal prism volume displacement indicate. The 
long-period effects at Station 3 might then be regarded as resulting from 
the longitudinal mixing of the waste discharge, whereas the similar effects 
at the other stations may be regarded as reflecting a combination of advec- 
live and diffusive processes. 

The diurnal components of the DO and BOD sj. »ctra may be interpreted 
as expressing the effects of diumal variations in waste discharge and In 
photosynthetic activity of the planktonic population. Diumal variations 
in waste discharge would affect stations within a tidal excursion of the 
outfall more strongly than those beyond this distance, whereas photosyn- 
thetic activity in the system would be more pronounced at those stations 
exhibiting the higher nutrient concentrations, generally reflected in higher 
mean stream BOD's. Interpreted on this basis, the diumal spectral com- 
ponents suggest that the extent of the upstream tidal excursion is between 
Stations 3 and 4, whereas Stations 6, 7, and 8 are subject to considerable 
photosynthetic activity in addition to the diumal waste load variations. 
The diumal components of the DO spectra at each station correspond in 
size to the respective BOD components; this is a result that might be 
expected from the theory of DO-BOD relationships in streams. 
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The semidiurnal component of these spectra reflects the advective mo- 
tion of DO and stream waste load due to tidal action. The magnitude of 
this component at each rtation is a measure of the longitudinal concentra- 
tion gradient of each parameter within a tidal excursion of the station. The 
DO spectra show large semidiurnal effects at Stations 3, 4, and 5, and 
relatively smaller ones at Stations 6, 7, and 8; these results indicate that 
there is a large DO gradient in the upper part of the system and a small 
one in the lower part Examination of the mean DO values at each station 
suggests the existence of an oxygen-sag regime in which Stations 3 and 4 
represent a zone of rapid degradation and Stations 5, 6, 7, and 8 a zone 
of critical DO and the beginning of recovery, a situation that is in close 
agreement with the spectral results. From this point of view the region of 
critical DO is the location that has the smallest semidiurnal spectral com- 
ponent, in this case Station 7. 

The semidiurnal BOD spectral components present a picture somewhat 
different from the DO spectra. At the two upstream stations this component 
of the BOD spectrum is small. There is a significant increase in this com- 
ponent at Station 5 and again at Station 6, whereas the semidiurnal effects 
at Stations 7 and 8 are of the same magnitude as that at Station 6. These 
results suggest a low BOD gradient in the upper reaches of the system and 
a high gradient in the lower reaches, and Station 5 represents a region of 
transition. The strong gradients near and below the waste outfall suggest 
that the semidiurnal spectral components are affected by variations in the 
initial mixing of the waste load throughout the volume of water passing 
the outfall in a tidal excursion. It is apparent that the mean BOD’s do not 
show such changes in gradient along the estuary. Comparison of th; BOD 
gradient regime shown by the spectral analysis with that shown by the 
mean BOD values at each station illustrates the sensitivity of spectral 
analysis as a tool in estuarine engineering. 

It has been the purpose of this discussion to demonstrate how spectral 
results can be interpreted in terms of familiar sanitary engineering con- 
cepts. It is not intended to suggest that the spectra can supply informa- 
tion in addition to that discussed. Spectral analysis is purely and simply a 
tool for the analysis of time-series data. It provides by itself no theoreti- 
cal insight into natural processes, but it does permit one to examine indi- 
vidual periodic components of the data, with a minimum of interference 
from other components. As with any other statistical technique, the final 
interpretation of these results must be based on an understanding of the 
natural process, not on some magic numbers produced by the manipula- 
tion of data. 
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DESIGN CRITERIA FOR SPECTRAL ANALYSIS 

I n the design of any spectral analysis program the requirements for pre- 
cision in each spectral estimate, for resolution of sufficient spectral bands, 
and for eliminating aliasing, or folding, must all be considered and bal- 
anced against each other before the sampling program and data analysis 
are begun. 

An acceptable balance among these requirements can be established by 
careful choice of the sampling interval, the record length, and the number 
of lags used in computation. 

The choice of appropriate sampling and computational factors should 
be based on reasonable assumptions of what the dominant periodicities in 
the system are. The shortest period it is necessary to resolve determines 
the sampling interval required, whereas the longest period necessary deter- 
mines the total record length. In most sanitary engineering applications, it 
may be assumed that diurnal fluctuations in waste discharges and In pho- 
tosynthetic activity will be of considerable importance. In tidal systems 
there will be significant semidiurnal periodicities; there may also be some 
effects of this short a period in waste discharge and in the biological 
regime for a nontidal system. As a general basis for experimental design, it 
may be assumed that resolution of 24-hour periods will be required in a 
nontidal system and that resolution of 12-hour periods will be required 
in a tidal system. 

The shortest period it is theoretically possible to resolve with a given 
sampling interval is the period that is twice the sampling interval. That 
Is, with a sampling Interval of 6 hours, it is theoretically possible to resolve 
a 12-hour period. From a practical standpoint it is not possible to do this, 
since the 1 12- hour estimate would be the shortest period computed, and 
there would arise the question whether this is a valid estimate of a 12-hour 
period or whether it merely represents the aliasing of periods shorter than 
12 hours. It is advisable that the sampling Interval chosen be small enough 
Ujprovuk spectral estimate* for several periods shorter than the expected 

“ ,hM "* p '" od h <***» •**«* •» "v 

* ru,e . 0 ^ *humb, it is suggested that a maximum sampling interval 
to teso,vc 24-hour periods and an interval of 4 
hours to resolve 12-hour periods. In general, a sampling interval of not 
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more than one-third the length of the shortest significant period is recom- 
mended. 

There are no clear criteria for determining the longest period that can 
be resolved from a given record. The longest period resolved in any par- 
ticular analysis (other than the “infinite” period estimate) is determined 
by the number of lags used in computation and by the sampling interval. 
In general, computation to a number of lags greater than 10 percent of the 
total number of measurements in the record is not recommended, i.e., for 
a record of 140 measurements, computation to no more than 14 lags is 
recommended. Each additional lag used in computation reduces the pre- 
cision of all spectral estimates computed, and is generally regarded that 
the 10-percent value affords an optimum balance between precision of indi- 
vidual estimates and resolution of spectral components. As a guide in 
determining the required record length for design purposes, it may be 
assumed that a record length at least 10 times as long as the longest signifi- 
cant period to be resolved will be required. For many sanitary engineer- 
ing field surveys, the 24-hour period is about the longest significant period 
it is necessary to resolve from field survey data; in such cases a minimum 
record length of 240 hours would be required. 

As a basis for selecting the number of lags to be used in computation 
10 percent of the number of measurements in the record is used as an 
upper limit. The number of lags finally chosen will probably be based on 
the resolution believed to be required. For example, if it is desired to 
separate diurnal and semidiurnal components, the number of lags chosen 
must be large enough to include several estimates between the 24- and 
12-hour estimates, so that there is a negligible amount of interference 
between the major components. An acceptable degree of precision is then 
obtained by increasing the record length if the numbers of degrees of 
freedom on which each estimate is based give a confidence band that is 
too large to give usable results. 

Any type of time-series record can be subjected to spectral analysis if 
it represents sampling at uniform time intervals and if there are no 
missing points; if a few points are missing, however, a limited amount of 
interpolation may be done. Interpolation of the mean value of all the 
measurements or linear interpolation of a missing point between two 
measurements is the usual approach. There are no definite criteria that 
serve as an indication of how much interpolation can be done in any 
particular case; a general consideration of the process suggests, however, 
that if the missing points are widely scattered, up to possibly 5 percent of 
the data may be interpolated without serious effects on the computed 
spectra. 

Data suitable for spectral analysis are obtained most simply and 
cheaply by means of automatic sampling and recording equipment. Con- 
versely, spectral analysis offers the most effective means of analyzing and 
correlaHng the large quantities of data produced by such instruments. 
This does not mean, cf course, that the computation of spectra is limited 
to time-series data obtained from automatic devices. The DO and BOD 
data discussed here were obtained by conventional manual sampling pro- 
cedures, whereas the tide height record was obtained from an automatic 
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recording device. If the requirements of uniform sampling interval and 
record length are met, the means by which the data are obtained is imma- 
terial. 

The computation of spectra is most cheaply and efficiently accom- 
plished by high-speed digital computer. The spectra presented here were 
computed on an IBM 704; the total cost, including preparation of the 
data for the computer, was estimated at less than 10 dollars per spectrum. 
Computation by hand for all but very short record lengths is prohibitive 
from the standpoint of time, money, and accuracy. 




INTRODUCTION TO CROSS-SPECTRA 

T he discussions in the preceding chapters were concerned entirely with 
the generalized harmonic analysis of a single data record. The end 
product of such an analysis is a spectrum of magnitudes of variance, showing 
the frequencies and magnitudes of the dominant periodic components of the 
record. 

Now consider the spectral analysis of a pair of records analyzed to- 
gether. The computations are quite similar, but the use of two records 
instead of one permits the calculation of some additional statistics, which 
have some quite direct and highly significant engineering uses. The means 
by which these statistics may be calculated it shown in some detail later; 
at this time let us briefly outline the types of information expected from 
cross-spectrum calculations and indicate what these may mean as aids in 
determining engineering parameters of estuarine systems. 

In the spectral analysis of a single record all phase relationshh * 
lost. For example, a record of dissolved oxygen (DO) might be pi 
by several different effects, which vary over a diurnal period. The 
change in water temperature might have an effect on DO, and the 
in photosynthetic activity of algae may also have an effect on D r 
of these effects would be diurnal in frequency, but the peak of DO 
photosynthesis could occur at a different time from that reficc 
change in water temperature. With but a single record to work fr 
the total magnitude of the resulting diurnal variation in DO can l 
mined. The phase relationships between these two diurnal effect 
be established from a single record. 

One of the powerful consequences of the cross-spectrum calcu 
the maintenance of phase relationships between dominant period 
ponents of the two records being analyzed. Cross-spectral am 
simultaneous water temperature and DO records would show V 
many hours the daily peak in DO follows the peak in daily ten 
and the size of the variation. (Since these two variables are ii 
related, the phase lag would show that they are out of phase, 
phase lag would be greater than 180 degrees.) Similarly, cross 
analysis of simultaneous sunlight and DO oxygen records would s 
the calculation of phase lag relationships between sunlight inter 
photosynthetic oxygen production. 



The use of two records simultaneously results in the computation of 
four sets of spectra: the individual power spectrum of each of the- two 
records, and two spectra resulting from the two cross-correlations of the 
two records, These two spectra, designated the cospectrum and quadrature 
spectrum t are components of the cross-spectrum itself. (The relationship 
of these two components to each other and to the cross-spectrum will be 
discussed in detail,) 

For each of the frequency bands for which spectral density estimates 
are made, there therefore result four estimates of spectral density reflect- 
ing the relationships of two parameters, From these four estimates several 
other extremely useful statistics can be calculated, 

Coherence is analogous to the square of a correlation coefficient. 

Phase lag is the angular time lapse between respective maxima or 
minima in the two records. 

Response junction is a quantitative estimate of the amount of variation 
in one record that is associated with similar variation in the other 
record. 

These statistics are extremely easy to interpret in terms of the physical 
environment through standard forms of mathematical analysis. 

Coherence has many of the properties of the square of a correlation 
coefficient and is often interpreted as such. In magnitude it is expressed 
on a 0 to 1.0 scale and is dimensionless, The coherence measures the 
closeness with which changes in one record follow those in the other. This 
closeness may be in both magnitude of variation and the closeness with 
which the frequencies correspond, If the one record is affected by several 
driving forces with the same dominant periodicities, the coherence may 
reflect that fraction of the variation in one record which is associated with 
the other, In concept, then, coherence tells us how closely the two records 
are related. 

Fhase lag is a measure of the time needed for changes in one record to 
bring about corresponding changes in the other. It is computed in terms of 
radians but can easily be converted to hours, days, or any other con- 
venient time unit. So, to go with the coherence, there is a measure of time 
relationships between two parameters, 

The response function is the spectrum of one record computed from the 
spectrum of the other record and from the cospectrum and quadrature 
spectrum of the two records. It is expressed in units of the original record 
and shows the total magnitude of the variance in the one record that is 
associated with the other. 

These three statistics, then, define three important factors in relating 
two parameters to each other. First, a measure of how closely associated 
they are: the coherence; second, a measure of their dynamic (or time- 
wise) relationship to each other: the phase lag; and third, a measure of the 
magnitude of the association: the response function. 

Examination of the relationship of one set of data to another is almost 
always concerned with establishing some sort of cause-and-effect relation. 
This is true in cross-spectral analysis as well as other types of analysis. 
Accordingly, one of the records to be subjected to cross-spectral analysis 
is termed the input and the other is called the output . 
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Figure 12, Typical Input-Output System 



The use of these terms is to regard the environment as a “black box.” 
That is, it is not possible to state precisely what is happening in the envi- 
ronment; but it is possible to measure what goes in and what comes out, 
and from these measures deduce a litile about what is happening in the 
enigmatic “black box” (Figure 12). In computing cross-spectra, then, it 
is assumed measurements are made of one causative factor and one 
resulting condition, and that one record is then an input and the other 
one an output. 

This terminology is a convenience in interpreting the results and a 
proper choice of input and output records can usually be made on physi- 
cal grounds. For example, in an analysis of a temperature record and a 
DO record, the temperature is the input and the DO is the output. 

The mathematics themselves are not affected by this choice, but the 
proper choice of input and output will avoid the trauma of trying to 
visualize negative phase lags, for example. 
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COMPUTATION OF CROSS-SPECTRA 

T he general scheme for the calculations of cross-spectra is the same as 
that for the computation of individual power spectra. There are some 
differences in computation and additional complexity of notation brought 
about by the use of two records instead of one. 

From an operational standpoint, the calculation of cross-spectra is 
much more involved than the computation of individual power spectra 
because the cross-spectrum scheme of analysis not only includes the com- 
putation of four spectra in each analysis but also includes further opera- 
tions with the resulting spectra. For this reason no attempt is made to 
follow through an actual arithmetic example to its conclusion. Instead a 
very simple algebraic example is used to illustrate the differences in com- 
putation and approach which are required in this scheme of analysis. 

Suppose that there are two records of data to subject to spectral ara’y- 
sis. In order to avoid having to carry a mean value throughout the discus- 
sion assume that the arithmetic mean is zero for each record. Let 

A — a record of 100 temperature observations at uniform time inter- 
vals, and 

B = a record of 100 dissolved oxygen observations of the same sys- 
tem over the same uniform time intervals. 

The individual sequential values of each record may be identified as 

Qi 9 & 2 i 03 , • • • flioo, and bit ^ 2 * ^3, • • • ^ioo* 

As the initial step in the analysis, form the autocorrelation of each rec- 
ord and the two cross-correlations of the pair of records. The autocorrela- 
tions for each record at lags 0 and 1 will follow this form: 



C 0 (A) = . . 


• +0ioo # 0ioo) YOO* 


(8) 


Ct(A) = (afa i +a 2 *a i +a 3 ra i + . . 




(9) 


Co(B) = (bi*bi-\-b 2 *b 2 -^- bs*bi-\- . 


• . +&ioo*&ioo) Jqq* 


(10) 


Ci(fl) = (.bi*b2-\-b2*b3-\-b3*bt-\- . . 


• . +^09*^ioo) £2 • 


(11) 
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From the pair of records the cross-correlations are formed. Here enters 
a complexity. Look briefly at how the cross-correlations are calculated. At 
lag 0, it is quite straightforward; the sequential values of each record are 
multiplied together, so that there is something of the form 

When the first lagged product is computed, however, there are two 
possible sequences of products which may result: 

( 1 ) 

(2) b\ % a2‘\-b2*ciz-\-bz*ai-]- . . . -f-Z>9o*tfioo 

It is immediately apparent that the sums of these two sequences of 
products will be the same only when the two records have identical values. 
When there are two different records, then, there are two cross-correla- 
tions. The cross-correlation functions used in the cross-spectrum scheme 
of data analysis actually consist of the sums and differences of the two 
basic cross-correlations. 

The general equations for these two cross-correlation functions are 



l K-r 

S, + = 2(n—r) £ (x,y, +r +x l+r y t ), 


(12) 


t n— r 




Sr = 2(n—r) £ ( x >yt+r-x t+ ,y,), 


(13) 



where x and y denote values of the records being used, 

t = 1,2,3, . . . = sequential index of values, 

r = 0,1,2, . . . ,/w = lag number, 
m = total number of lags. 

For the example, at lag 0 there are the equations 
$o + = 2(M) U^x # ^i+«i^i) + (^2+^2)+ ... 

+ (tfl00*frl00 + d f 100*f>10o) ] 

= 2(100) — a — a2*by) + ... 

+ (^ 100*^100 — tfl00 * frl0o )] 

For lag 1 there are the equations 

Sl H — 2 ? 99\[( a l*^2+ fl 2*^l)+ (^2 # ^ 3 +^2)+ . . • 

$1“ — ^2*^1 ) +^3*^2)— ' 0362)+ • . . 

+ (' 00 — 0100 *^ 99 )] 

The next step in the analysis is to compute the Fourier transforms for 
the two autocorrelations and for the two cross-correlation functions. For 
the cross-correlation function S f +, the Fourier cosine transform (C r ) is 
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computed using the relationships 



i Ti m—1 -i 

Co = — [j (So + +S m +)+ 13 S ' + J’ 

0 f 1 m — 1 rri 'l 

Cr = — | y [$>++(— l) r d’m+H- £ 5,+COS j, 



1 ( . ’ m — 1 ^ 

= i \ i tS 0 +-f-<— £ (-1)%+ J. 



(14) 

(15) 

(16) 



For the cross-correlation function S r , the Fourier sine transform 
Qr is computed according to the relationships 



Qo = 0, 


(17) 


0 m — 1 




q; = - £ s t - sia*z, 
m q % m 


(18) 


£ 

II 

© 


(19) 



This step performs the same smoothing function for the cross-correla- 
tion functions as the Fourier transformation does for the autocorrelation. 
The use of the cosine transform on one of the pair and the sine transform 
on the other is a key step in the calculation which permits the computation 
of phase lags from the results of the analysis. 

The reason for this lies in the fundamental geometrical relationship of 
the sine and cosine. That is, in a given right triangle, the sine of one of the 
angles is the ratio of the opposite side to the hypotenuse, while the cosine 
is the ratio of the adjacent side to the hypotenuse. Numerical values of 
the sine and cosine are equal for angles which are 90 degrees out-of-phase 
with each other, so that the effect of multiplying one cross-correlation by 
the sine and the other by the cosine is to weight each cross-correlation by 
a function which suppresses components of the variance which are not 90 
degrees out-of-phase with each other. 

The final step in the treatment of the cross-correlation functions is a 
weighting process exactly analogous to that used in the computation of 
individual power spectra. The results consist of the cospectrum (C r ) and 
the quadrature spectrum ( Q r )> which are computed using these relation- 
ships: 



0.54 (Co'+C,'), 


(20) 


0.23 Cr—i -f-0.54 Cf-['0,23Cr+i , 


(21) 


0.54 


(22) 


0.54 ( Qo+Qn , 


(23) 


0.23 0r — l ”l-0.54 0r'-fO.23 Qr + l\ 


(24) 


0.54 (Q m ~ i )• 


(25) 



These two results consist of two spectra which are 90 degrees out of 
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phase with each other. They may be considered components of the output 
spectrum. A crude geometrical analogy to what has been done is presented 
in Figure 13. Instead of a resolution of vectors into components, there is 
a resolution of spectra into components. 

The four basic statistics computed in this analysis are these: 

Uri = po./er spectrum of input record 
U ro = power spectrum of output record 
C , = cospectrum 
Qr = quadrature spectrum 

A number of other statistics can be calculated from these basic results. 
There are considered here three of these statistics which appear to have a 
basic relationship to studies of the dynamics and water quality of estuarine 
systems. 



OUTPUT 




INPUT 



Figure 13. Geometrical Analogy to Cross- spectral Analysis 



The coherence (H r ) is computed according to the relation 

C , a +( 21 



H r = 



( 26 ) 



Hr is a dimensionless number, and is generally regarded as analogous to 
the square of a linear correlation coefficient. The terms in this equation 
are directly related to terms in standard formulations for a linear correla- 
tion coefficient. One such form is 




J , 2 

V 







= linear correlation coefficient of y or x 



= covariance of y and x values 
= variance of x values 
= variance of y values 
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Hr represents a component-by-component expression for such a statistic, 
while r ty represents an averaged value for the entire records ot x and y 
values. 

The phase lag (£ r )is computed by the relation 

L, = arc tan (2?) 

The units are radians, which can be converted into any convenient time 
units related to the sampling interval. Since Q, is a sine transform and 
C r is a cosine transform, the reason for this relationship is obvious. 

The response function spectrum (R r ) is computed by the relation 

R r = H r >U ro (28) 

and has the units of the output spectrum. 

These three statistics in addition to the four spectral values form the 
basis for a coherent scheme of data analysis to interpret sampling results 
in terms of the physical environment. 




INTERPRETATION OF CROSS-SPECTRA 

T he purpose in subjecting data to any form of mathematical manipula- 
tion is to establish relationships between parameters of the physical envi- 
ronment. Any statistic is useful only when it can be interpreted in terms 
of the environmental situation. It is necessary, therefore, to relate the 
statistics computed in cross-spectral analysis to properties of the environ- 
mental system. 

The interpretation of individual power spectra has already been dis- 
cussed in some detail and needs no further comment here. 

The cospectrum and quadrature spectrum are the basis for several fur- 
ther calculations showing the degree of correlation between two records. 
The actual values of the cospectrum and quadrature spectrum are right- 
angle components cf the response of the output spectrum to the input 
spectrum. In terms of the physical environment, these components have 
no direct significance; their value and utility lie in the fact that from them 
can be calculated coherence, phase lags, and response spectra. 

It is worth noting that the relative magnitudes of the cospectrum and 
quadrature spectrum govern the phase lag, and their absolute magnitudes 
dominate the coherence. 

Coherence . The coherence is roughly analogous to the square of a cor- 
relation coefficient and, in the simplest use of spectral analysis, it is often 
interpreted exactly as a linear regression correlation coefficient is inter- 
preted. Because of the means whereby it is calculated and because it is 
derived from environmental rather than isolated and controlled measure- 
ment, it may be interpreted in a number of other ways peculiar to environ- 
mental studies and to cross- spectral analysis. 

Conceptually, the term “coherence” is precisely descriptive of what the 
statistic means. Coherence is a measure of how coherent the two records 
are; that is, how similar the variances of the two records are at each 
spectral estimate and how well the fluctuations in the one record follow 
those in the other. The coherence, therefore, is not really merely a measure 
of linear correlation of two variables; it has implications concerning both 
the magnitudes of variances observed and the shape of the spectral com- 
ponents examined. 
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This point may be clarified with an example. Consider two records for 
cross-spectral analysis: a record of sunlight intensity and a record of dis- 
solved oxygen determination. The actual data records for several clear 
days might look like Figure 14. 

The solar intensity record reflects the fact that, when the sun sets, the 
paiameter immediately reaches its minimum value, while the DO record 
reflects the fact that there is a sequence of events and processes between 
the sun’s energy and its ultimate effect on DO. 

These two records are quite dissimilar in shape, not necessarily because 
of a lack of correlation between the two, but because of the nature of the 
processes involved. The coherence calculation reflects only the actual 
numbers in the two records, so the coherence computed between these 
two records might be a bit lower than the physical situation actually war- 
rants. As with any other statistic, coherence must be interpreted with a 
clear realization of how the number was derived. 

As a general rule of thumb in interpreting coherencies, it can be stated 
that a coherence is always lower than, or equal to, the true degree of 
correlation between two records. 

One of the most useful ways to interpret coherence is to regard it as 
the fraction of the variance in the output record that is related to similar 
variance in the input record. For example, if the coherence between the 
diurnal components of the solar intensity and DO records was 0.60, we 
could state that about 60 percent of the diurnal variance in DO is ex- 
plained by the diurnal component of variance in solar intensity. 

This interpretation of coherence follows directly from its means of 
computation. The coherence, as calculated, is the ratio of the sum of 
the cross-covariances to the product of the individual power spectra; that 
is, the coherence is actually a ratio of squared variances. Its interpreta- 
tion as a fraction of variance is therefore quite logical. 

Response function spectrum. This statistic is the calculated response 
spectrum of one record to another. The response function spectrum shows 
what the output spectrum would be like if the input record were the only 
parameter dominating it. The response function spectrum has the same 
physical meaning as the individual power spectra have, and may be 
interpreted in much the same manner. 

The response function spectra of one output to several different input 
parameters can be computed separately and the results compared to deter- 
mine what major driving forces are affecting any particular parameter. 
Response function spectra do reflect the phase relationships between the 
input and output records; in this respect they may therefore be inter- 
preted more advantageously than the power spectra of the individual 
records. 

Phase lags. The phase lag between two records estimates the time delay 
in the process or processes which the two records measure. The physical 
significance of this delay may be interpreted in several ways depending on 
the nature of the records and processes. 

. For records of different parameters at a single station the phase lag may 
quantitatively estimate the time required for a biological, chemical, or 
physical process to occur. For example, a phase lag between air tempera- 



42 



ture and water temperature would estimate the time required for the 
transfer of thermal energy from the air to the water or vice versa ; the 
phase lag between solar intensity and DO would measure the time delay 
in the biological and chemical processes of photosynthetic oxygen pro- 
duction. 

For records of the same parameter at different stations, the phase lag 
may be a measure of advective transport between two stations. For ex- 
ample, the phase lag of tidal height between two stations might show the 
progression of the tidal wave up an estuary. The phase lag between DO's 
at two stations might, however, measure the change in DO due to con- 
sumption or reaeration during the time of the phase lag as well as advec- 
tive transport; a considerable degree of caution in interpreting phase lags 
in this manner is indicated when a combination of physical, chemical, and 
biological processes is possible. 

In any form of interpretation of cross-spectral results the procedures of 
the calculation must be clearly understood and their significance kept 
firmly in mind, 

The cospectrum and quadrature spectrum are variances calculated from 
the mean values of the records, just as the individual power spectra are 
calculated. This does not affect the values of coherence, phase lag. or 
response spectrum, but it can affect the interpretation that may be placed 
upon these values. For example, if the mean DO of a record is very near 
the value of DO saturation for the temperatures and salinities concerned, 
one would expect the autocovariance function and the cross-covariance 
functions to be smaller than they would be if the mean DO of the record 
was about 50 percent of saturation, because the DO saturation value is a 
limiting value for the record. The same is true if the mean value of the 
record is very near 0 percent of saturation. 

The relationship of the mean value and limiting values of the recorded 
parameter to the record variance and other descriptive statistics may be 
illustrated with some data collected on the Potomac estuary. 

In Figure 15 are shown mean values and variances of DO records with 
145 sequential values at nine stations near Washington, D. C. This is the 
fresh-water part of the estuary and the mean temperature at each station 
was between 28.0 and 29.0 degrees centigrade. The mean DO results 
therefore represent comparable values. While the actual shape of the curve 
sketched in Figure 15 may be debated, the relationship between mean DO 
and magnitude of record variance is unmistakable. 

Low record variance may be associated with both unpolluted and highly 
polluted reaches of a river or estuary, white high DO variances would be 
associated with the midrange of mean values where the processes of de- 
oxygenation and reaeration are competing and about equally effective. The 
frequency distributions of the records would be skewed toward the mid- 
range of values if the mean values are near zero or saturation, since these 
represent limiting DO values, while the distribution would probably be 
Gaussian for means near 50 percent of saturation. 

The spectral estimates from records near these limiting values would 
also be biased, even though the calculation of the spectral statistics would 
not show it. It is apparent, then, that a scheme of interpretation in which 
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the spectral estimates were assumed to rcsutt from a sinusoidal fluctuation 
about the mean (or any other type of uniform fluctuation) could lead to 
quite erroneous results. 

There is one more statistic which may be computed from cross- spectral 
results and which is peculiarly useful in environmental studies. This Is the 
response, which is actually a rate statement of the process between two 
records. Response is a very useful but somewhat dangerous statistic to use, 
because its dimensions and name imply that a great deal more is known 
about a process than may actually be the case. 






Response is calculated from the response function spectrum and the 
input spectrum by dividing the spectral estimates of the former by those of 
the latter and then taking the square root. Another form of estimation of 
the response is to divide the sum of the estimates of the response function 
spectrum by the similar sum of the input spectrum and then take the 
square root, 'this latter procedure is more logical from an engineering 
standpoint since what is actually observed in nature is a mixture of the 
dominant periodic component of both records; this procedure is equiva- 
lent to saying that the standard deviations of the response function spec- 
trum and the input spectrum are valid means of quantifying the response 
of the one record to the other. 

This approach is valid providing the input and output records are statis- 
tically similar. For example, the frequency distributions of both records 
should be alike: either both Gaussian, or both skewed in the same direc- 
tion, or both of similar kurtosis (peakedness). If the two records are not 
statistical^ similar, then the use of a response calculation based on them 
may lead to deceptive results. 

The response has units of the output record divided by the input record: 
ppm DO/ppm BOD, ppm DO/cfs river discharge, ppm BOD/lb sewage 
load. These rates are, of course, based on variance from the mean values 
of the input and output records. While they may actually represent funda- 
mental rates for the system being studied, their interpretation should be 
based on the observed mean values unless there is other contributing infor- 
mation available to support a different interpretation. 
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SPECTRAL ANALYSIS AS A DATA 
PROCESSING SYSTEM 

T he major purpose in obtaining and analyzing any set of environmental 
data is to determine relationships between the parameters measured. 
This purpose presupposes some vague idea of how the system of interest is 
operating and what the controlling parameters are. Figure 16 is a gen- 
eralized flow chart of the dissolved oxygen system in an estuary which 
shows conceptually how DO may be expected to respond to environ- 
mental changes. 

Such a general conceptual appreciation of the workings of a system is 
a prerequisite for any coherent scheme of data processing and analysis. 
This does not mean that the internal complexities of » ach part of the sys- 
tem must be known; it does mean that the major interrelationships of the 
system should be understood. 

One of the great advantages of spectral analysis as a scheme of data 
analysis is that extremely useful quantitative results can be obtained with 
only the most general knowledge of how the driving forces and the water 
quality parameters of interest are related. 

For example, a record of solar radiation can be analyzed with the DO 
record for the same period, -and information on the effect of sunlight on 
DO can be obtained without direct knowledge of the photosynthetic oxy- 
gen production chain. The results obtained are, of course, purely empiri- 
cal numbers based on analyses of field data; they do not necessarily rep- 
resent any new theoretical insight into the system. They are, however, 
answers about relationships in a particular system at a particular time and, 
as such, are extremely useful in developing an engineering appraisal of an 
estuarine water quality problem. 

Spectra] analysis is most effective and useful when the environmental 
system under study is dominated by highly variable and complex driving 
forces. In estuarine systems the complexities introduced by the harmonic 
nature of tidal and diurnal driving forces make spectral analysis especially 
useful. 

The scheme of spectral analysis as outlined here represents a direct 
coherent system for the analysis of environmental data. Where records of 
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data from automatic sampling devices are available, routine processing of 
the records by spectral analysis would provide rearly all the statistical 
evaluation of these data that would ever be needed on a routine basis. 

With such a scheme of data processing and analysis available, it should 
be possible to analyze in depth on a routine basis the records from auto- 
matic sampling stations. This would reduce by a very large fraction the 
number of man-hours spent in unprofitable data analysis. A routine data 
analysis procedure in which basic environmental records, such as tem- 
perature, solar radiation, river dischaxge, and tidal height, were run against 
water quality parameters such as DO would prove a continuous evalu- 
ation of the response of any estuarine system to the environment and 
establish quantitatively what environmental factors are significant in any 
particular system. 

This is not to say that spectral analysis should be used in all estuarine 
situations. In cases where manual sampling must be used to measure 
certain parameters, the situation must be carefully evaluated to deter- 
mine if the expense of manually collecting data for spectral analysis is 
justified because the problem at hand is too complex for any other means 
of analysis to work. This is a matter for individual evaluation. 

No scheme of data analysis, particularly one based primarily on statis- 
tical manipulation of data, can substitute for engineering knowledge and a 
logical approach to an environmental problem. Spectral analysis is a 
scheme of data processing which can produce valid statistical results 
rapidly and in highly complex situations; the relationship of these statis- 
tics to the physical environment must, however, be determined through 
engineering logic. 

This chapter contains five examples of the use of spectral analysis in 
analyzing and interpreting time series data in river and estuary field 
studies. The first concerns the tracing of a conservative pollutant in a 
river, the second concerns the dilution and decay of a nonconservative 
pollutant in a river, the third illustrates the estimation of diffusion of a 
nonconservative pollutant in an unslratified tidal river (/.*., negligible 
salinity), the fourth shows how river flow affects temperature and OO in 
a tidal system, and the fifth deals with the relationship of river flow to 
stratification in a stratified estuary. 

The design of any spectral analysis computation requires the balancing 
of several conflicting computational factors; the rationale for choice of 
design criteria is presented in detail for the first of these examples to 
illustrate differences in governing criteria for design, and to demonstrate 
the flexibility of spectral analysis as a statistical tool in extracting differ- 
ent types of information from time series data. 

Example 1. Tracing a Conservative Pollutant in a River 

A discharge of calcium chloride in the wastes from a large chemical 
plant was affecting the raw water hardness of a municipal water supply 
some 200 miles downstream. Between these two points was a large 
reservoir and two large tributary inflows, both of which also had Urge 
impoundments. 



Table l.a. Input Record. Salt Waste Discharge 
(CaClg to River 1000 lbs/ day) 



Year 


Jan. 


Feb. 


Mar. 


Apr 


May 


June 


July 


Aug. 


Sepl. 


Oct. 


Nov. 


Dec. 


1950 


3496 


4137 


3067 


2023 


4781 


1938 


2603 


1951 


2745 


1584 


1690 


4725 


1951 


2686 


2478 


2866 


2376 


2217 


1834 


1737 


1789 


2143 


1416 


1858 


4228 


1952 


3518 


1975 


2548 


1877 


1632 


1203 


972 


993 


888 


971 


2139 


1668 


1953 


1874 


2255 


2946 


1861 


3596 


1765 


2431 


3260 


1013 


809 


898 


1786 


1954 


6469 


4759 


2484 


1970 


1934 


1283 


3457 


1359 


1261 


1282 


1218 


2084 


1955 


1480 


2653 


2604 


1930 


1478 


1325 


1209 


1339 


1049 


963 


1306 


1360 


1956 


1460 


1385 


1892 


2725 


2123 


1672 


1528 


1301 


1239 


1267 


1157 


2229 


1957 


2712 


3152 


1835 


1659 


1349 


1165 


1074 


808 


1272 


929 


2304 


1999 


1958 


1693 


1661 


1461 


1695 


1782 


1463 


818 


1300 


1015 


936 


1172 


1769 


1959 


1788 


1298 


1984 


2165 


2115 


1288 


1193 


1556 


1 05 1 


2250 


2979 


2048 



Information derived from direct sampling of stream water quality, raw 
water supply hardness, and waste discharges showed that, at a mean flow 
of 10,200 cubic feet per second over a four-year period, 26.2 percent of 
the total hardness in the municipal water supply was due to this waste 
discharge, and from sampling during a subsequent one-year period with 
a mean flow of 13,600 cubic feet per second, 20 percent of the hardness 
was due to this waste discharge. 

Confirmation of these results was desirable, and cross-spectral analysis 
of the waste discharge quantities with raw water hardness was used to 
obtain this confirmation. 

The data available consisted of the two time series records showm in 
Tables l.a and l.b. The values in each table are monthly averages of 
the parameters measured over the same consecutive 10-year period. 

Observe that one record (Table l.a) is in units of thousand pounds/ 
day of calcium chloride while the other (Table 1 .b) is in units of mg/1 of 
hardness. The dimensions carried by records used for cross-spectral analy- 
sis need not be identical for computational purposes. Some of the resulting 
statistics will, however, carry the combined units of the original records; 
attention should be given, before calculation, to the most useful units for 
expression of the final results and any necessary conversion made imme- 
diately. 

In the present example a response in terms of mg/1 hardness per lb/ 
day CaG* would result; this is a useful number for the purposes of the 
analysis, so no conversion is needed and the records can be used as they 
arc. The record length and sampling interval arc already established by the 
data available, so that the only design criterion which can be altered is 
the number of lags used in computation. Since there are 1 20 points in each 
record, the rough 10 percent rule-of-thumb suggests that 12 lags can 
be used readily. Recent practice in spectral analysis suggests that a much 
larger number of lags can be used (as much as one-third of the record 
length has been used by the author, even on very short records, with good 
results). 
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Table l.b. Output Record. Water Supply Hardness (mg/l) 



Year 


Jan. 


Feb. 


Mar. 


Apr. 


May 


June 


July 


Aug. 


Sept. 


Oct. 


Nov. 


Dec. 


1950 


102 


89 


75 


74 


103 


104 


106 


102 


90 


102 


180 


123 


1951 


132 


115 


72 


91 


79 


105 


87 


104 


126 


141 


130 


129 


1952 


95 


93 


73 


76 


103 


116 


109 


124 


125 


125 


134 


138 


1953 


118 


90 


104 


102 


74 


m 


126 


127 


128 


141 


132 


139 


1954 


87 


77 


64 


6$ 


62 


88 


122 


114 


139 


128 


137 


120 


1955 


120 


95 


79 


86 


86 


88 


97 


118 


116 


119 


120 


120 


1956 


130 


70 


97 


87 


94 


106 


96 


114 


125 


132 


141 


109 


1957 


117 


76 


78 


80 


78 


81 


79 


90 


92 


93 


88 


86 


1958 


94 


80 


75 


75 


81 


93 


93 


108 


108 


109 


115 


131 


1959 


95 


72 


110 


93 


84 


66 


94 


87 


95 


88 


106 


101 



To decide the number of lags to use in this example, these factors were 
considered: 

1. What will be the major periodic fluctuations in the system? 

Since the waste discharges are being carried along in the river flow, 
the major periodic components in hardness variation at *he water 
supply will be closely related to the hydrologic cycle. This cycle has 
a very strong annua] component, and because of the large it tpound- 
ments on the river and its major tributaries, much of the shorter 
period fluctuations in flow will be damped out and the annual cycle 
in flow will govern the downstream hardness variations. It is de- 
sirable then, to choose a number of lags in which the annual cycle 
will be shown clearly. 

2. How much resolution is needed? 

High resolution, which is obtained by using a large number of lags, 
is needed where there may be periodic components of nearly the 
same frequency in the data. For example, if we suspected that there 
were significant variations in the data *hich recurred at 10-month 
intervals as well as at 12-month intervals, we would go to a high 
enough number of lags so that these two components would be sep- 
arated in the analysis (this point will be illustrated subsequently). 
In the present case the existence of the strong annual cycle and the 
damping effect of the dams suggests that the annual cycle will domi- 
nate the entire spectrum, and that other periodic fluctuations in the 
flow and In the hardness will be suppressed. 

If there is doubt about what periodic components may be present 
in any set of data, the best approach is to compute power spectra 
on the individual records at several different lags, and from these 
results observe what periodic components are significant; the final 
design of analysis can then be made with considerable confidence. 

3. How much statistical reliability is needed? 

The higher the number of lags used, the wider is the confidence 
band for each estimate. When a purely random process governs a 
system, this is a very significant problem; in a system which is 



highly deterministic the statistical measures of confidence limits are 
much less important. Even so, it is good practice to keep the number 
of lags used in the analysis as low as is compatible with the resolu- 
tion needed. 

The actual choice of what number of lags to use is determined by 
selecting several possible numbers of lags, calculating the recurrence 
periods for which these choices will give spectral estimates, if necessary 
calculating the degrees of freedom of each estimate and the confidence 
interval, and then making a choice among them. 

For the present case 6, 12, 15, and 18 lags all seemed to be equally 
reasonable. The equation for calculating the recurrence period corre- 
sponding to each lag number was given earlier (Equation 6) and is 

Tr = ( 6 ) 

where 

T f = period corresponding to the rth lag, 

M = time interval between data points (I month here), 
m = total number of lags, 
r = lag number, o < r < nu 

For each of the choices shown above, the recurrence periods at which 
spectral estimates will be obtained were calculated and are presented in 
Table 2. 

The 6-lag, 12-lag, and 18-lag choices all will provide spectral estimates 
for a 12-month period. The 15-lag choice will give spectral estimates for 
15 month and 10-month recurrence periods, but none for 12 months. 
Th.s means that a strong 12-month spectral component would be averaged 
between the 1 5-month and 10-month estimates if this number of lags is 
used; this is not desirable, so this choice is rejected. 

The 6-lag choice has a 12-month period estimate, but it is adjacent to 
the long-period (*) estimate. If either record has a variance component 
greater than one year in recurrence, then the estimate of this component 
will overlap the 12-monlh estimate and it will be impossible to separate 
them. Therefore the 6-lag choice is ruled out. 

The 12-lag choice has one estimate between the long-period and 12- 
month estimates, and the 1 8-lag choice has two. The higher resolution of 
the 18-lag choice is at the cost of some loss of statistical confidence in the 
result. 

The numbers of degrees of freedom for each spectral estimate of the 
12-lag choice and the 18-lag choice were calculated using the method 
given by Blackman and Tukey, which is based on the Chi-square test. This 
calculation proceeds in this manner: 

(1) Calculate effective record length ( Tm ) from actual record 
length (7 ) and total number of lags used (m): 

r „ ^ r - 1/3 (m) = 120 - 1/3 (12) = 116 

(2) Calculate equivalent band width (IF r ) for each estimate 
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Table 2. Recurrence Periods at Various Choices of Total Lags 



Lag Number 


Period in Months at Each Lag Number for Total Lags of: 




6 


12 


15 


18 


0 


00 


oe 


00 


00 


1 


12.0 


24.0 


30.00 


36.00 


2 


6.0 


12.0 


15.00 


18.00 


3 


4.0 


$.0 


10.00 


12.00 


4 


3.0 


6.0 ' 


7.25 


9.00 


5 


2.4 


4.8 


6.00 


7.20 


6 


2.0 


4.0 


5.00 


6.00 


7 




3.4 


4.30 


5.10 


8 




3.0 


3.75 


4.50 


9 




2.7 


3.33 


4.00 


10 




2.4 


3.00 


3.60 


it 




22 


2.70 


3.30 


12 




2.0 


2.50 


3.00 


13 






2.30 


2.80 


14 






2.10 


2.60 


15 






2.00 


2.40 


t6 








2 25 


17 








2.10 


IS 








2.00 



(3) Calculate change in frequency between spectral estimates (A/) 

_ J ! 

J 2r„ (2) (116) 

(4) Calculate number of degrees of freedom (k) for each estimate 

k _ _ (2) (U6) _ 19 , 

This is very roughly equivalent to saying that a person could have as 
much confidence in the magnitudes of each spectral estimate as he could 
have in the average value of 20 pieces of data. For the 18-lag computa- 
tion a similar calculation gave 12.3 degrees of freedom, or a confidence 
level equivalent to that of having an average based on 13 pieces of data. 
Since both lag choices have the resolution required, and the numbers of 
degrees of freedom are both within reasonable working ranges for water 
pollution engineering, the choice between the 'wo becomes a subjective 
one. 

The 12-lag choice was made, primarily because compulation time 
would be considerably shortened w hile a sufficiently high degree of reso- 
lution was maintained. 

The cross-spectral analysis computations were then made using the salt 
waste discharge record as the input and the water supply hardness at the 
output. Table 3 is a representation of the computer output obtained by 
cross-spectral analysis erf these two records. The computer program was 
designed so that the output contained not only the individual power spectra 
and the cross-spectrum of the two records, but also many subsidiary cal- 



5 ) 



Table 3. Computer Output of Cross-Spectral Analysis: 
Salt Waste — Water Hardness Problem 

CROSS-SPECTRAL ANALYSIS 

INPUT RECORD — SALT WASTE — 1000 LBS/DAY CACL2 
OUTPUT RECORD — WATER SUPPLY— MG/L WATER HARDNESS 
120 VALUES 12 LAGS 1.0 TIME INTERVAL 



INDIVIDUAL SPECTRA CROSS-SPECTRA COMPONENTS 



LAO 


INPUT 


OUTPUT 


COSPECT. 


QUADSPECT. 


0 


L744E 03 


6.322E 01 


5.294E 02 


0 . 


i 


1.622E 03 


I.013E 02 


-9.280E 02 


-I.287E 03 


2 


J.336E 03 


I.322E 02 


-2.234E 03 


-2.746B 03 


3 


1.198E 03 


6.999E 01 


--9.244E 02 


-1.616E 03 


4 


7.130E 04 


I.383E 01 


2.196E 02 


-3.057E 02 


3 


5.054E 04 


1.193E 01 


6.219E 01 


-3.762E 01 


6 


4.434E 04 


1.192E 01 


-2.477E 01 


1.71 IE 02 


7 


4.677E 04 


LOOSE 01 


2.202B 01 


1.282B 01 


8 


4.009E 04 . 


6.863B 00 


— 3.3S4E 01 


— 1 .29 1 E 02 


9 


3.177E 04 


7.180E 00 


-7.806E 01 


-8.924E 01 


10 


2.730E 04 


8.875E 00 


-2 802E 01 


— 1.1 50E 02 


II 


2.916E 04 


9.839B 00 


— 5.617E 0! 


— I.010E 02 


12 


3.267E 04 


I.026E 01 


— I.314E 02 


0 . 


LAO 


PERIOD 


COHERENCE 


RESPONSE 


PHASE (DEO.) 


0 


LONO 


0.0234 


1.606E 00 


0 


1 


24.00 


0.1531 


1.551E 01 


234.2 


2 


12.00 


0.6175 


8.163E 01 


230.9 


3 


8.00 


0.4135 


2.894E 01 


240.2 


4 


600 


0.1436 


I.986E 00 


305.7 


3 


4.80 


9.0087 


1.040E 01 


328.8 


6 


4.00 


0.056 6 


6.75 IB 01 


98.2 


7 


3.43 


0.0014 


L402E 02 


30.2 


8 


3.00 


0.0647 


4.440E 01 


255.3 


9 


2.67 


0.0616 


4.424E 01 


228.8 


10 


2.40 


0.0579 


5.147E 01 


256.3 


11 


2.18 


0.0464 


4.57IE 01 


240.9 


12 


2.00 


0.0515 


5.280E 01 


180.0 


LAO 


COHERENCE ROOT 






0 


.139 








1 


.391 








2 


.786 








3 


.642 








4 


.379 








3 


.093 








6 


.238 








7 


.037 








8 


.234 








9 


.248 








10 


.241 








11 


.215 








12 


.227 









OVERALL RESPONSE = .01162 



culations which ate very useful, but are time consuming to calculate with- 
out the use of the computer. Some of these are not needed in the present 
example, and will be explained when included in the other examples pre- 
sented here. 
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Figure 17 show! the individual power spectra of each record, as well as 
the response spectrum of the pair of records as also presented in Table 3. 

The waste discharge spectrum in Figure 17.a shows a pronounced long- 
period component and a “probably" strong annual component. The wt>rd 
"probably" is used because there could be a strong 24-month component 
and the annual component would reflect the overlapping estimates of the 
annual and 24-month components. However, the variance estimates at 
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8-month and 6-month periods drop rapidly and consistently from the 
magnitude of the 12-month estimate, so it is probable that the magnitude 
of the 24-month spectral estimate is the result of overlap between the 
long-period and 12-month estimates. If what was desired from anlysis of 
this record was the spectrum of the waste discharge, this analysis would 
be inadequate and there would be two courses of action open: 

First, repeat the analysis at a higher number of lags to ascertain the 
magnitude of the annual component by computing estimates at several 
periods between the long-period and annual periods. (See f he preceding 
discission on the higher resolution obtained by using 18 lags instead of 
12 bgs.) 

Second, detrend the record by running a linear (or second-order) re- 
gression of waste discharge against time, compute deviations from the 
regression line instead of the mean of the entire record, and then compute 
the spectrum on these deviations. This is a potentially misleading proce- 
dure, since the investigator is committed to the assumption that the regres- 
sion line actually is a deterministic input into the record, instead of a 
merely statistical device for making something out of nothing. 

Fortunately, in the present example, the actual spectrum of the waste 
discharge is important only as it relates to the hardness of the water supply 
downstream. 

The spectrum of the water supply hardness, Figure 17.b, shows a very 
pronounced annual component, a smaller amount of long-period variance, 
and very small variance recurring at other periods. It is apparent that the 
12 lags in this case have been quite satisfactory in resolving and quantify- 
ing the variance of the water supply hardness record. The same approach 
used in the further analysis of the waste discharge record would be appli- 
cable for this record if additional information were required from the 
analysis of this one record. 

Figure 17.c shows the response spectrum of the water supply hardness 
record to the waste discharge record. The dominant feature of this spec- 
trum (which results from a calculation using the cross-spectra and coher- 
ence of the two records) is the annual component of variance; all other 
variance components are of negligible magnitude. Such a strong response 
at one recurrence period and a lack of response at other periods suggests 
that the assumptions about the importance of the impoundments, in 
smoothing out all except the annual flow variations, were quite correct; 
these results also suggest that the water supply hardness is a highly deter- 
ministic function of the salt waste discharge. 

In Tabic 3 are shown the square roots of the coherence for each 
spectral estimate. The annual recurrence period has a coherence root, 
roughly equivalent to a correlation coefficient, of 0.8, which is a very high 
result for this type of analysis.* The coherence reflects the amount of 

• While no detailed analysis of the relationship between coherence root and cor- 
relation coefficient has been made as yet, the theory of probability suggests that 
coherence root is analogous to the correlation coefficient for a Joint probability 
density function, and it would therefore be the product of two correlation coeffi- 
cients; one relating variarxe magnitude, and one relating consistency of periods. If 
this analogy fs accepted, then the coherence root of O S would he equivalent to the 
product of two such correlation coefficients, each of about 0.9. 
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correlation between the magnitude of variance in input and output records 
for each recurrence period, as well as the consistency in the repeating of 
the period of recurrence itself. In this case, if there were some variation 
in the time of annual peak runoff from year to year, there would be a 
decrease in coherence even if the magnitude of hardness variation were 
the same. 

The results presented in Table 3 contain the information necessary to 
calculate the amount of hardness in the water supply attributable to the 
salt waste discharge. 

The response spectrum is in the same units as the output spectrum, (i.e., 
ppm hardness) 2 ; this spectrum is what the output record spectrum would 
be if the input record atone generated the output. The single estimate unit 
response , i.e., ppm hardness per 1000 Ib/day salt discharge, is the square 
root of the ratio of each response spectrum value to each input spectrum 
value for each of the recurrence period estimates. The overall unit response 
is the square root of the ratio of the summed response spectrum values to 
the summed input spectrum values and has the same units as the single 
estimate unit response. 

Both of these are valid statistics, each particularly informative in cer- 
tain types of situations. 

The single estimate unit response, R s , is most useful when there is no 
transfer of energy from one recurrence period in the input to another 
recurrence period in the output. For example, if the major recurrence 
periods in both the salt discharge and water supply hardness records in 
the present example had been the annual cycle, then the assumption log- 
ically would be made that changes in the annual cycle in the input would 
bring about changes in the annual cycle of the output record, and 
would be a suitable response statistic to use. 

The overall unit response, R 0 , is used when there is evidence that there 
has been a transfer of energy from one recurrence period in the input spec- 
trum to a different recurrence period in the output spectrum. In the pres- 
ent example there has been a transfer of energy from the long-period 
estimates in the spectrum of the input record to the annual estimate in the 
spectrum of the output record. This transfer of energy between periods is 
the result of reservoir storage of salt discharge for many months within 
the river system, followed by its release from storage when the annual 
hydrologic cycle brings about a partial flushing of the stored salt dis- 
charges, which is a major controlling factor in downstream water hard- 
ness. For the present example, then, R 0 is the appropriate response statis- 
tic to use. 

This statistic is given directly in Table 3, the computer printout, as 
overall response; for :he present example, R 0 = 0.0116 ppm hardness/ 
1000 lbs. per day of salt. 

The average percent of salt hardness at the water plant directly attribu- 
table to the salt discharge may be calculated directly from R 0 in this 
manner: 

j ^ * R 0 X mean salt discharge 

mean salt hardness (%) — - — • (29) 

mean hardness 

In the somewhat oversimplified “mathematical moder of salt hardness 
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given above, there has been made the implicit assumption that the river 
discharge has no effect on the water plant hardness. The salt hardness 
estimates using stream water quality data were, however, based on the 
observed condition that total hardness varied inversely with river dis- 
charge and on the conclusion, from this and other data, that tributary 
inflow of salt was in direct proportion to river flow and had little effect on 
total hardness concentration. 

The total hardness is, therefore, affected by the varying dilution of the 
salt waste in the river, an inverse function of river flow. 

The cross-spectrum calculation gives results which are averaged for the 
entire period of record and are related directly to the mean values of all 
controlling conditions during the period of record. In this case R 0 = 0.012 
is the overall unit response at a mean salt discharge of 1,935,000 lbs/day 
and a mean water plant hardness of 102 ppm. It also applies for the mean 
river flow of 12,450 cfs during the period of record. To include the dilut- 
ing effect of the river flow on the salt waste, the earlier expression for 
salt hardness is adjusted by including a dilution factor term in the clas- 
sical form: 

VJ5 m = KS nf (30) 



where V m and V„ are river discharges and S m and S n are corresponding 
salt hardness concentrations. 

This can be included in equation (29) to give the result: 



where 



S* 



( Vm ) _ RW ( V m ) 

(V n ) ~ H (Kn)’ 



(31) 



W = salt discharge for survey n (lbs/day), 

H = mean hardness (ppm), 

V m = mean river discharge (cfs), 

S m ~ mean salt hardness concentration (ppm), 

V„ = river discharge for survey n , 

S n — salt hardness for survey /?. 

Table 4 summarizes comparable calculations of salt hardness for each 
of the two field surveys and for the mean of the two surveys. 

From the stream water quality data obtained at two different periods, 
corresponding salt hardness percentages of 26.2 and 20 percent were 
obtained, an average of 23.1 percent which compares well with the 22.8 
percent calculated with equation (31) in this manner: 

For the entire period of record . . . 

_ (.012) (1935 ) (12450) _ „ R 
102 (12450) 

For the two surveys . . . 

* = (.012) (1852) (12450) _ _ 

106(10240) 

, = (.012) ( 1722) (12450) = - 

98(13600) y 
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AH of the differences between corresponding salt hardness values shown 
in Table 4 are within the limits of precision of the chemical determination 
for hardness. 

The cross-spectrum calculation thus provides an excellent check of the 
other calculations. The spectral analysis is based only on the salt waste 
discharge and the water plant hardness, not on the stream survey data, so 
the two methods of calculation are independent. 

The cross-spectral analysis of these two records provides another sta- 
tistic not obtainable by the other means of calculation used in the analy- 
sis of these data. This is the mean time of passage of the salt waste from 
the point of discharge to the water plant. 

Table 3 shows the phase lags associated with each recurrence period. 
The annual recurrence period is a major feature in the spectrum of the 
output record and also has the highest coherence between records; this 
annual cycle reflects the yearly flushing of the impoundments during high 
spring river flows. The 7.7-month phase lag shown in Table 3 for this 
annual cycle is an average time of passage of salt waste between the two 
points of record; in this case this number probably represents reservoir 
residence time to a large extent. 

Knowledge of time of passage can be used as a pollution control meas- 
ure through timing of releases to minimize effects at downstream points. 

This example shows how information in available plant records can be 
extracted through use of spectral analysis. In the present case not only 
did the analysis of plant records alone provide the same information, 
which, using conventional techniques, required a stream survey in addi- 
tion to the plant records, but it also provided additional useful informa- 
tion. 

A slightly more refined approach to this particular problem would 
include the cross-spectral analysis of river discharge and water plant 
hardness records to compute an overall unit response factor for hardness 
to river discharge; this factor would be a more precise measure of the 
relationship of hardness to river discharge than the dilution effect only, 
as assumed here. 

Example 2. Dilution and Decay of a Nonconservative Pollutant In a 
River 

The first example showed how the effects of a conservative pollutant 
can be examined quantitatively using spectral analysis. In the present 
example, a similar approach is used in the evaluation of the effects of 
two sewage discharges on a river. This example is based on the standard 
sanitary engineering concepts of the relationship of DO and BOD, with 
the additional use of the sunlight-DO photosynthesis relationship exhi- 
bited in Figure 16 and quantified through the cross-spectrum calculations. 

In this example, the starting point is the waste load of the sewage treat- 
ment plants as described by their BOD loads; the impact of these loads on 
the stream is traced through the changes in stream BOD and DO as they 
vary with waste loads, with sunlight, and with river flows. 

The physical situation is that shown in Figure 18. The two waste out- 
falls were sampled at three-hour intervals to obtain 127 consecutive 
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10.25 miles) outfall a 



STATION I tO O MILES) 

STATION 2 (1.25 MILES) 
OOTFALL B 11.70 MILES! 

STATION 4 12.75 MILES) 



STATION 6 15.25 MILES) 



STATION 8 (10 35 MILES! 



Fiyurt 18. Schtmatic Diagram Showing Loc«l»n« of OutfaMa and Sampling Station* 




l samples (about 16 days). Figures 19 and 20 present graphically the data 

| obtained in this sampling program. Observe that the waste discharge from 

i Outfall A is about 10 times as great as that from Outfall B. The spectra 
of the two waste discharges are shown in Figures 21 and 22. 
s Comparison of volume and concentration spectra with the total load 

spectra from each of these two outfalls shows that the total load of BOD 
! discharges is not closely related to the pattern of volume or concentration, 
f This is of interest because of the tendency to base estimates of the flue- 
f 
I 

I 
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Figurt 21. Spectra of W**»« Oitchvgt Parametert, Outfall A. 




tuation in a sewage plant wasteload on the fluctuations in volume. Total 
loads are used in the further computations here. 

The total load spec'rum of Outfall A exhibits long-period and diurnal 
components of about the same magnitude, while there is also a very small 
semidiurnal component. If it is assumed that all the periodic fluctuations 
are sinusoidal (the dangers of this assumption hive been discussed), the 
amplitude of each of these components can be estimated from the spec- 
tral estimates.* The sine-wave amplitudes (SWA) corresponding to these 



* The variance of a sine wave from its mean value is equal to times the 
square of its amplitude. ^2 



64 



(lb B0D/DAY)2 x 10* (ppm B0D)2 (mgd)? 






Figure 22. Spectra of Wtsta Oiacharge Pararrmn, Ou tf»ll B, 



\ variances are: long-period, 58,000 lbs. BOD fi /day; diurnal, 53,000 lbs. 

BOD 6 /day; semidiurnal, 16,500 lbs. BOD 6 /day. The mean was 136,000 
lbs. BODs/day. The three sets of components indicated account for about 
91 percent of the total variance of the record. 

The total BOD 6 load spectrum of Outfall B exhibits a large long-period 
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component, SWA 4700 lbs./day, and semidiurnal components, SWA 
5500 lbs./day. These account for about 79 percent of the total record 
variance. The mean waste discharge is 15,100 lbs./day. 

Both discharges exhibit wide fluctuations in total load, and the fluctua- 
tions In the Outfall B load are proportionately greater than those in the 
Outfall A load. Outfall A discharges through a diffuser across the river 
channel, while Outfall B is a point source on one shore. These two facts 
are of importance in discussing the character of the stream BOD responses 
to these two outfalls. 

Figure 23 presents the spectra of stream BOD and the total load re- 
sponses of stream BOD for both outfalls at five stream stations. Station 1, 
which is above both outfalls, exhibits primarily long-period components, 
and only quite small shorter-period effects. From Station 2 downstream 
the effects of the strong diurnal components in the waste discharges are 
quite apparent. 

Several environmental conditions, in addition to the spectrum of the 
entering waste load, determine the BOD response spectrum at any stream 
station. 

First, the relative magnitudes of the BOD load already in the stream 
and the BOD of the entering waste load can determine whether the BOD 
response of the stream is large enough to be observed and what is the 
magnitude of the overall response. 

Second, the turbulent-mixing regime and the initial dispersion of the 
waste in the stream exert some control over the response spectrum. Tur- 
bulent mixing in the stream tends to eliminate small variations in concen- 
tration first and larger variations later; ultimately responses from all sizes 
of variations become smoothed out into secular trends in the data. An 
intense turbulent-mixing regime would tend to bring about this transition 
more rapidly than a weak one. 

Third, the nonconservative nature of the BOD parameter means that 
there is a time-dependent decreasing trend in the BOD response at stream 
stations. 

Fourth, there may be some indirect response of the stream environment 
to the waste load, and this may exert some effect on the stream BOD 
response. For example, nutrients entering stream in the waste load may 
stimulate alga! growth; the increased algal population is then measured as 
part of the stream BOD at downstream stations. 

In the present situation, the discharge from Outfall A represents about 
22 percent of the total load in the stream below the outfall, whereas the 
discharge from Outfall B represents about 2.4 percent of the total load 
in the stream below that outfall. The response spectra to the two dis- 
charges are, however, very similar in magnitude at Station 4 as well as 
at Stations 6 and 8, In Figure 22 the diurnal response components to 
Outfall A are also about the same at Station 2 as at Station 4. 

The lack of change in the magnitudes of the diurnal components of the 
response spectra indicates that the longitudinal turbulent-mixing regime 
in this stretch of the river is not strong enough to affect the initial disper- 
sion of wastes from this outfall. There is a similarity in magnitude of spec- 
tral components in response to Outfalls A and B at Station 4. A well- 
mixed input from Outfall A and a poorly mixed input from Outfall B 
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affects the stream BOD at Station 4, and it is logical that the relative 
magnitudes of the response spectra at this station should bear little resem- 
blance to the relative magnitudes of the waste discharges. 

The change in mean value brought about by the waste inflow, any 
secular trends in the waste discharge, the sampling and determination 
errors in the BOD test, and the longitudinal turbulent-mixing regime gov- 
ern the long-period response components. Turbulent-mixing phenomena 
lend to cause a shift in the dominant response components toward the 
longer periods and ultimately these effects appear as a secular trend in 
the mean value at a particular sampling station. When a nonconservative 
parameter, such as BOD, is the criterion considered, the biochemical 
changes occurring at the time may obscure the shift in the spectral re- 
sponse. 

The long-period response components at Station 2 probably reflect 
secular trends in the discharge for Outfall A, whereas the smaller ones 
at Station 4 reflect the decreased importance of the long-period fluctua- 
tions in waste discharge from both outfalls in the presence of the higher 
mean stream BOD at this station. 

The response spectra at Station 6 show significant changes from those 
at Station 4. Suppression of the diurnal response components to both 
outfalls suggests that the controlling phenomenon is the turbulent mixing 
between Stations 4 and 6, whereas the increase in magnitude of the long- 
period response components to Outfall B indicates that there has been 
some shift in response to longer periods as the waste discharge has be- 
come mixed into the stream. These may be related to the difference in 
initial dispersion of the two waste discharges, since the response to Cut- 
fall A shows no similar effect. 

At vStation 8, the diurnal response components to both outfalls are of 
about the same magnitude as those computed at Station 6, indicating 
either a somewhat-less-intense turbulent-mixing regime in thL section of 
the river than was indicated in the reach between Stations 4 and 6 or a 
limit to the effectiveness of turbulent-mixing processes as far as longitu- 
dinal waste dispersion is concerned. An interesting phenomenon in the 
response spectrum to Outfall A is the appearance of a semidiurnal re- 
sponse component. This component does not appear at Stations 4 and 6, 
and its appearance so far downstream suggests that it may be an indirect 
result of the waste discharge rather than a direct one. Inorganic nutrients 
introduced in the waste discharge, that tend to increase the algal popu- 
lation, may explain this component. If this effect does exist, it may affect 
the entire spectrum at this station 

The rationale presented in the preceding paragraphs serves as the 
physical basis for a quantitative interpretation of the results of the cross- 
spectral analysis of these data. Tt>e total response, f.c., the sum of all 
the periodic response components, decreases as stations farther from each 
outfall are reached. This decrease is a function of the longitudinal tur- 
bulent-mixing processes, which tend to smooth out concentration fluctua- 
tions resulting from the waste discharge, and of the rate of biochemical 
oxidation of the total waste load in the stream. It is not possible with 
present knowledge to separate these two effects In the response spectra 
to obtain independent estimates of biochemical degradation and turbu- 



Table 5. Correlation Coefficients of Response Spectra 
to Outfall Discharges 



Recurrence 



Correlation coefficient of response spectrum 



spectral 

component, 

days 


2 


to Outfall A at Station 
4 6 


8 


to Outfall B at Station 
4 6 8 


X 


.61 


.37 


.40 


.16 


.48 


.72 


.18 


3 


.55 


.38 


.37 


.15 


.55 


.59 


.12 


1.5 


.66 


.73 


.49 


.52 


.43 


.28 


.23 


1.0 


.81 


.88 


.60 


.71 


.80 


.57 


.62 


0.75 


.80 


.82 


.48 


.64 


.61 


.31 


.57 


0.60 


.67 


.42 


.13 


.44 


.13 


.35 


.25 


0.50 


.66 


.33 


.39 


.80 


.19 


.46 


.12 


0.43 


.54 


.19 


.41 


.71 


.21 


.40 


.11 


0.38 


.36 


.14 


.39 


.35 


.30 


.25 


.19 


0.33 


.37 


.15 


.37 


.43 


.42 


.16 


.36 


0.30 


.16 


.23 


.58 


.21 


.48 


.12 


.29 


0.27 


.13 


.18 


.38 


.26 


.49 


.08 


.22 


0.25 


.05 


.40 


.15 


.30 


.69 


.10 


.27 



Correlation coefficient s V coherence. 



tent mixing from the spectral results alone. The spectral results do, how- 
ever, provide phase lag relationships that can be interpreted as time-of- 
passage information, and biochemical degradation rate constants have 
been determined from tong-period BOD’s on this system. From this in- 
formation the decrease in response due to biochemical loss in the waste 
load can be determined, and then some estimate can be made about tur- 
bulent mixing and dispersion in this reach of the river. 

Phase lags can be calculated for each periodic component of the re- 
sponse spectrum as shown in Example I. The best phase lag to use for 
interpretation as a time of passage is obtained from the response com- 
ponent that has the highest degree of correlation with the waste discharge 
spectrum. Table 5 shows the correlation coefficients computed for each 
spectral response component at each station for both outfalls. While these 
correlation coefficients cover a wide range of values, those relating the 
diurnal responses are fairly high, above 0.8 for Stations 2 and 4. The 
somewhat lower correlation coefficients at Stations 6 and 8 might be 
expected since mixing of the waste into the river, which may be charac- 
terized as a random process, leads to distortion and ultimate disappear- 
ance of the short-period fluctuations related to the waste discharge. 

Table 6 presents the phase lag in hours of the diurnal components of 
stream BOD, as well as responses calculated from spectral results. 

The phase lags show the time elapsed between passage of the diumal 
peak past the outfall sampling point and its passage past the stream sam- 
pling station. Time of passage in the river itself may be obtained by 
subtracting phase lags to each successive sampling station. Tie negative 
phase lag indicated for Outfall B at Station 4 is an artifact of the compu- 
tation resulting from the fact that the diumal peak in BOD from Outfall 
B does not coincide with the diumal peak in BOD at Station 4; it has no 
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physical significance in the interpretation of phase lags as times of pas* 
age in the river. Table 7 presents time of passage calculated from these 
phase lags. The agreement between times of passage calculated from each 
outfall is fairly good. Velocities for this reach of the river in these dis- 
charge ranges have been estimated at about 1.0 miles per hour from 
volume displacement calculations; this is also a good check of the spec- 
tral results. 

The deoxygenation rate constant for BOD in this reach of the river is 
assumed to be 0.16 day“ A at 20 # C for the purpose of estimating the loss 
in BOD due to biochemical oxidation. At the mean temperature of the 
water during this study, the deoxygenation rate constant would be about 
0.24 day" 1 . From this rate constant and the time-of-passage estimates of 
the preceding example, the total response at each station in the absence 
of any loss through turbulent mixing can be calculated. A measure of the 
turbulent mixing can then be obtained from the remaining change in total 
response at each station. The pertinent statistics incorporating this 
approach are in Table 8. 

The statistic of particular interest is the dilution factor, D, defined as 
the ratio of total response (when the BOD loss is due to biochemical loss 
only) to the observed total response. If it is assumed that each sampling 
point is representative of a volume of the water passing that point, then 
the value of D is regarded as a volume-of-dispersion ratio. For example, 
if a unit quantity of the waste discharge from Outfall A is dispersed in a 
volume V as it passes Station 2, then, as this same quantity of waste passes 
Station 4 it is dispersed in a volume equal to 1.18V. The similarity of 
values of D for Outfall A at Stations 4 and 6 to those for Outfall B at 
Stations 6 and 8 is very striking. This similarity suggests that the turbulent 
diffusion characteristics of waste dispersal in this river are more closely 
related to the concentration of the waste at specific points rather than to 
some intrinsic property of the receiving stream. 



Table 8. Dispersion of Wastes Between River Stations 



Station 


Waste from Outfall A 


Wade from Outfall B 




Number 


Rt 


Rl 


D 


Rt 


Rl 


D 


2 


.733 


_ „ 


__ 


_ ^ 







4 


.589 


.693 


1.18 


.538 


— 


— 


6 


.379 


.557 


1.47 


.424 


.500 


MS 


8 


.326 


.348 


1.07 


.264 


.264 


1.47 



Rt ~ total response (BOD mg/1 ) 

Rl “ total computed response for BOD loss due to decay only 



D = dilution factor r= — 



Although this form of interpretation and calculation of turbulent diffu- 
sion processes does not have the intrinsic beauty found in calculations of 
eddy diffusion coefficients, it does have the virtue of providing an insight 
into how the waste material disperses in the stream, thereby giving an 
appreciation of how the waste assimilation capacity of a stream can be 



utilized most effectively to dispose of waste effluents. It is worthwhile 
noting that, as shown in Table 8, turbulent diffusion is more significant in 
reducing the BOD observed at successive stations than are the natural 
assimilation processes, which are the usual bases for estimating a stream’s 
waste assimilation capacity. 

Spectral calculations offer a fresh approach to the problems of turbu- 
lent diffusion and waste dispersion in streams. This brief discussion aims 
to show* only how' these problems yield to the techniques of spectral analy- 
sis. No attempt is made to explore these possibilities fully. 

Example 3. Diffusion of a Nonconservative Pollutant in an Unstratified 
Tidal River 

The two preceding examples showed the use of cross-spectral analysis 
in rivers and streams. In these cases the dominant harmonic components 
W'ere the diurnal variations in sewage discharges and in the photosynthetic 
production of oxygen. 

The present example introduces a major harmonic component common 
to all estuarine systems — the motion associated with the tide. The situa- 
tion presented involves sewage discharges entering an unstratified tidal 
river, so that the effect of tidal movement can be examined free of the 
complications of salinity-dominated stratification. 

This example is an extension of the analysis of the Potomac Estuary 
near Washington used in Chapter 3 to illustrate the interpretation of indi- 
vidual power spectra. The present discussion is based on the cross-spectra 
of these BOD and DO records and the use of the solar radiation record 
in conjunction with them. 

The spectrum of solar radiation obtained during the field study is 
shown as Figure 24. The strong diurnal component is expected, but it is 
interesting to note the appearance of a significant semidiurnal component. 

The source of this component is not known, but several possibilities 
exist. These data were taken at a Weather Bureau research station and 
precautions were taken to avoid back-radiation from local night lights or 
other local environmental artifacts. The semidiurnal component may be a 
moonlight effect or reflection of night lights of the city from clouds. The 
means of analysis would regard a situation in which there was a daytime 
peak of radiation (sunlight) and a nighttime peak of radiation (moonlight 
or city lights) as a semidiurnal effect of small magnitude. Whatever the 
cause, it still represents an impingement of energy on the instrument and 
consequently on the water. 

This discussion centers around the cross-spcctral analysis of three types 
of records: solar radiation at one point in the area, DO at nine river sta- 
tions and BOD at nine river stations. The cross-spectra computed with 
these records were: 



Input Output 

Solar Radiation BOD 

Solar Radiation DO 

BOD DO 



Solar radiation and the effects of waste load (including diffusion and 
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advection of wastes) are the primary driving forces affecting the DO in this 
system. While the stream BOD is a measure of some of the effects of the 
waste load, it is also a measure of other things. The algae population, 
which produces and consumes oxygen, also contributes organic matter to 
the stream BOD. The DO-solar radiation cross-spectra will show the net 
1 j effect of solar radiation on DO, which would include both the direct oxy- 

gen production or consumption and the change in BOD due to changes in 
algal cell material, f.e., decayable organic matter. 

The individual power spectra of the stream DO and stream BOD at the 
nine stations on the Potomac estuary are shown In Figures 25 and 26. A 
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schematic representation of the sampling station locations was presented 
in Figure 11. 

The DO response at each stream station to solar radiation is shown in 
Figure 27. The increasing importance proceeding downstream of diurnal 
fluctuations and longer variations is quite obvious. The suppression of 
semidiurnal response in comparison to the strong semidiurnal components 
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of the DO spectra themselves is also noteworthy. The diurnal variations j \ 

are due to differences in algal oxygen production and consumption in the 
daytime and at night, while the long-period effects result from the in- , 

creases and decreases in algal populations as they respond to changes in 
the emironment. 

The relationship of DO response to solar radiation with stream BOD 
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may be examined by plotting the total response, i.e., the sum of the re- 
sponses for the individual spectral bands, at each station as a function 
of the mean stream BOD at each station. This has been done in Figure 
28. (The point representing Station O is anomalous because of the ex- 
tremely high silt turbidity of the Anascostia River.) 

The sigmoid curve sketched in Figure 28 is based on the assumption 
that the changes in DO response are the results of changes in algal popu- 



7 ? 






(\ 



— n»i . _ . . 

/ 

! 




f f ljuf» 21 R»lr)Oft#ilp ftt too M D O. ^Mp«AM to totvr fttdtottoA 

ERJC ’* 

i 



i 



lation, and that the algal population is related to and Is part of the stream 
BOD. This type of curve Is generally typical of the response of a dodu- 
lation to a favorable environment. Measurements of the algal population 
during this study indicated that the numbers of algae were higher at lo- 
tions with higher stream BOD and that the numbers of algae fluctuated 
diumally, Part of this relationship is, of course, an artifact of the BOD test. 

This graph has demonstrated that the stream BOD and the solar ra* 
diation response of DO are relateo, and that the solar radiation is affect* 
ing the stream BOD values. (It may be safely assumed that stream BOD 
does not affect solar radiation.) 

It would be helpful to discuss the solar radiation and stream BOD as 
separate driving forces on the stream DO so that quantitative estimates 
of the effects of waste discharges on DO can be made. To do this it is 
desirable to determine the extent of the the dependence of stream BOD on 
solar radiation and form some basis fo" eliminating this dependence 
from the total DO response to solar radiation or to BOD. 

The response of BOD to sunlight was computed from the stream BOD 
at each station and the solar radiation. Since there Is no reverse interaction 
of stream BOD on sunlight, we can assume that this response represents 
an additive effect as far as stream BOD variance is concerned, thus the 
effect of sunlight on stream BOD can be eliminated simply by subtracting 
the stream BOD response-to-sunlight spectrum from the stream BOD 
spectrum. This will provide a "corrected” stream BOD spectrum which 
can be regarded as a primary driver independent of sunlight. The net 
effect of solar radiation on DO is measured with the DO-solar radiation 
cross-spectrum; this includes the effect of the part of the stream BOD 
variance due to solar radiation. 

The use of the DO response derived directly from the stream DO- 
solar radiation cross-spectra has certain great advantages over attempt- 
ing to derive a DO response factor by using plankton counts, temperature, 
and stream BOD separately. For example, when the water temperature 
is highest the oxygen saturation value is lowest. This, however, might also 
be the time when algal oxygen production is highest, so that these two 
effects would tend to balance each other. While a partial correction for 
temperature could be made by using percent saturation of DO (at least in 
fresh water), this might actually obscure some relationships which could 
depend on the actual concentration of DO. It can be seen, therefore, that 
to follow the sunlight-DO system through the individual members of the 
chain would require a detailed knowledge of the amounts of interaction 
and the periods involved. While the direct computation of stream DO 
response to sunlight does not provide knowledge of the internal workings 
of this system, it does show the overall result; for many engineering appli- 
cations this is all that is necessary. 

The corrected stream BOD spectra are shown in Figure 29. These spec- 
tra are primarily the result of the quantities of wastes discharged into the 
estuary and how they mix into the receiving waters of the tidal river. The 
DO response spectra to corrected stream BOD can now be calculated by 
multiplying the DO response spectrum to BOD by (1— H r ) for the solar 
radiation-BOD cross-spectrum. This uses the coherence as the fraction 
of variance related to a particular driving force. In this case the assump- 
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tion is that all variance in stream BOD not related to sunlight is caused 
by waste discharge. The DO response spectra to corrected stream BOD 
are shown in Figure 30, 

These results may be interpreted in terms of the waste discharge and 
receiving waters regime in the Potomac Estuary. 

The waste volume discharges are all dominated by a strong diurnal 
fluctuation, while the tidal river motion is dominated almost entirely by 



80 



D.O. VARIANCE (ppm 2 x1(>4) 



1000 








760 








600 


— 




\ 


260 






\ 










0 


STATION 0 


STATION 1 


STATION 2 









ST A VI ON 6 STATION 7 STATION 8 

PERIOD OF RECURRENCE (HOURS) 

Figure 30. D.O. Response to Corrected Sveam BOO 

a semidiurnal rhythm. These conditions might be expected to have these 
effects on the stream BOD spectra: within the limits of a tidal excursion 
from an outfall, there may be a noticeable diurnal effect as well as a 
pronounced semidiurnal component. The relative magnitudes of these two 
components would depend upon the size relationships of the waste dis- 
charge, receiving water volume, and stream BOD already present and 
upon the intensity of mixing in the receiving medium. Outside the limits 
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of a tidal excursion from a major outfall, the stream BOD fluctuations 
would be dominated by the semidiurnal effect to the practical exclusion of 
diurnal effects. The magnitude of the semidiurnal effect observed at any 
station would depend on the longitudinal BOD gradient weighted by the 
tidal excursions; this effect might be regarded as a measure of the longi- 
tudinal turbulent mixing. 

The BOD spectra for Stations 3 through 8 exhibit dominant long-period 
and semidiurnal effects; while the diurnal effect is noticeable at Stations 
6, 8, and possibly 4, it is quite suppressed. From these spectra ' \ might be 
deduced that the volume of waste discharged is quite small compared to 
the* volume of the river into which it is discharged; the strength of the 
semidiurnal peaks would indicate the longitudinal BOD change is fairly 
small in the waters which flow past Stations 3, 4, and 5 during the tidal 
cycle, while the changes in the waters flowing past Stations 7 and 8 are 
rather large. Station 6 shows an intermediate condition. Figure 31 shows 
the volumes of water in the system at mean high water and at mean low 
water. The head of tidewater is assumed to be at Three Sisters Islands. 
This figure may be used to obtain a crude idea of the length of a tidal 
excursion in this reach of the system. For example, water opposite the 
Blue Plains outfall (Station 6) at high tide was about halfway between 
Stations 7 and 8 at low tide, while water at Station 6 at low tide lies 
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between Stations 4 and 5 at high tide. It should be clearly understood that 
| this type of estimate is quite crude and is based only on volume displace- 

ment. Not only are there large inherent errors In the volume measure- 
| ments, but also no account is taken of the current structure of the system, 

| the nature of the tidal wave (tranbiatoiy or standing), or of the effect of 

i river discharge. This figure is presented to give omy some crude idea of 

what the tidal excursion might be like. It can be seen from Figure 31 that 
\ the water flowing by Stations 3, 4, and 5 during a tidal period comes pii- 

1 manly from the region between Station 7 and a short distance above Sta- 

tion 3; some water would also come from the Anacostia River. The bulk 
of the water involved would be represented by the mean BOD values at 
i Stations 0, 4, 5, 6, and 7, which are 3,15, 3.34, 3.73, 4.10, and 4.30 

‘ respectively. Only a small amount of the water involved would be repre- 

) sented by the mean BOD of 2.44 at Station 3. It may be also noted that 

the BOD’s from Stations 6 through S would begin to show the effects of 
! larger tidal excursions, which would appear as a larger BOD change in 

j the water moving past the station. (If the longitudinal BOD changes in 

/ the tidal excursion at each station were weighted according to the length 

l of the tidal excursion, the results might provide some measure of the tur- 

bulent-mixing processes.) 

These results and the postulated mechanism may now be considered in 
/ comparison with the DO response spectra in Figure 30. It is immediately 

/ evident that the DO response (i.e., rate of change) to the diurnal BOD 

fluctuations is proportionately greater than to the semidiurnal fluctuations; 
) in some cases the diurnal effects are about the same magnitude as the 

j semidiurnal, as at Stations 6 and 8. Since the diurnal BOD component is 

! regarded as that part of the total stream BOD which has been most re- 

f cently discharged from the treatment plants, the large DO response to 

this component can be regarded as representative of the rapid biochemical 
\ utilization of waste materials well broken down by the treatment processes, 

while the gradual assimilation of the remaining materials is reflected in the 
j long-period and semidiurnal effects. 

f For estimation of waste assimilation capacity it is necessary to calcu- 

\ late the total DO response to stream BOD. This is done by dividing the 

; total DO response to BOD by the total corrected BOD variance at each 

j station; the square root of this quotient is the desired factor. For the data 

\ discussed here, these results were obtained: 
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What these factors say is that every ppm of BOD entering this section 
of the estuary will require about 0.5 ppm of DO. This factor is a measure 
of oxygen utilization by wastes in this part of the estuary and aJso a 
measure of flushing time for these wastes, if the deoxygenation rate con- 
stant is known. 

These calculations with similar analyses of river discharge, tide height, 
waste load, and temperature records as input with the DO and BOD rec- 
ords as output provided the basis upon which a detailed evaluation of the 
effect of river discharge on the waste assimilation capacity of the Potomac 
Estuary could be made. 

It was found that the waste assimilation capacity of the estuary in- 
creases with increasing river flow due to an increase in reaeration capa- 
bility combined with a decrease in estuarine temperature which increases 
DO saturation values in parts of the estuary where rapid waste assimila- 
tion is taking place (e.g., Stations 2 and 3). 

The details of these calculations are presented in the report * from 
which these results were taken. 

Example 4. Relationship of Change in River Flow to Temperature and 
Dissolved Oxygen in an Unstratified Tidal System 

The Potomac estuary is used again to illustrate how spectral analysis 
provides a mechanism for estimating the magnitudes of cause and effect 
relationships by indirect means. In this case, with river discharge the only 
record available for the inflowing stream, and with temperature and DO 
records on the estuary but not on the river, it was possible to determine 
the effect o f river flow on temperature and DO in the estuary, and even 
to estimate the DO concentration and temperature of the river above the 
estuary. 

Temperature affects the saturation concentration value for DO in water 
as well as the reaction rates of seme of the biological, chemical, and physi- 
cal processes into which DO enters. This discussion deals primarily with 
the relationship of temperature to the DO saturation value. In an estuary 
the water temperature will be governed by two major parameters: the 
advection of heat into the system in inflowing water, and the absorption 
of thermal energy by the water from the air and from solar radiation. The 
DO response to temperature will then result from the action of the two 
primary driving forces, solar radiation and river discharge, if this is the 
only major source of inflowing water. 

In Figure 32 the stream DO response function spectra to water tem- 
perature are plotted for each of the nine sampling stations. These spectra 
exhibit significant long-period, diurnal, and semidiurnal components at 
nearly all stations. It is useful to compare these spectra with those of the 
primary driving forces and of the water temperature itself. In Figure 33 
are presented the spectra of solar radiation, river discharge, and tem- 



* “Report on the Potomac River Basin Studies. Part VII, Technical Appendix.” 
U.S. Army Corps of Engineers, 1962. {Part VII was prepared by the USPHS, R. A. 
Taft Sanitary Engineering Center, Cincinnati, Ohio, and Region III, USPHS, Char- 
lottesville, Virginia.) 
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figure 331 Spectra erf Driving Form and Water Tamparatura 

peraoire for three typical stream stations, and temperature response spec- 
tra to river discharge for the same stations. 

The spectrum of solar radiation is primarily diurnal while that of river 
discharge is almost entirely long-period. The water temperature spectra 
are also entirely long-period, and all of the temperature response to river 
discharge is long-period. 

From this inspection of these spectra some conclusions can be drawn 
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about the DO response to temperature. First, it is highly unlikely that 
river discharge is directly related to anything but long-period response 
components; second, diurnal and semidiurnal components of water tem- 
perature are relatively unimportant; third, the diurnal components of the 
DO responses to water temperature are probably closely related to tidal 
flow and their magnitudes are more closely associated with the longitudi- 
nal gradients of temperature and DO than anything else. 

The magnitude of the diurnal component of DO response to tempera- 
ture is probably the result of a number of biological and chemical proc- 
esses related to decomposition of organic matter, algal activity, and oxy- 
gen saturation values. We do not know enough about the mechanisms 
involved to analyze these components further. The diurnal components of 
the temperature spectrum are so small that a more precise determination 
of the diurnal change in DO due to temperature could be obtained only 
with much longer records than are available here. 

The long-period components of the DO response spectrum lend them- 
selves to a somewhat more detailed analysis. River discharge is the con- 
trolling parameter for these components of the DO response spectrum, 
so that the statistic of interest now is the DO response to river discharge. 
The river discharge has two effects on DO in the estuary: (1) it warms 
or cools the estuary, thereby changing the potential oxygen saturation 
values in the estuary; (2) it advects into the estuary an amount of DO 
which is directly proportional to its flow rate and inversely proportional 
to its temperature. 

These effects may be examined quantitatively through cross-spectral 
analysis of the pertinent records. The volume of the river discharge is 
the only thing known about the fresh water inflow: neither the DO con- 
centration nor the temperature is known. It is interesting, therefore, to 
note that the effects of river discharge on the temperature and DO of the 
estuary can be determined without knowing anything about the tempera- 
ture and DO of the river discharge itself. 

Table 9 contains the water temperature spectra and the response spectra 
needed for further calculations. 

The phase lags are not included in any of these tables. It is noted that 
the temperature response to river discharge is out of phase; that is f de- 
crease in river discharge brings about an increase in temperature. This 
result would be expected for August, when these data were taken. At this 
time of year, water resident in the estuary would absorb more heat and 
lose less heat to the environment than would water tumbling down from 
the mountains in a free flowing stream. 

In Table 10 are presented some results calculated from the spectral 
results in Table 9. 

Columns 2, 4, 5, and 7 are the long-period spectral components (long 
+ 96 hour) taken from Tables 9.b, 9.a, 9.c, and 9.d respectively for each 
of the stream stations. 

Column 3 is obtained from Column 2 by dividing by the variance of 
river discharge (0.0373), and entering 0.25 times the square root in 
Column 3. (The factor 0.25 is a conversion from gage height to cfs for 
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the river discharge.) For Station 0 this computation is 



1 / 1.93 

4 \ 0.0373 



1.80 



Column 8 is obtained from Column 7 in the same way. 

The entries in Column 6 are the square roots of the quotients of Col- 
umn 5 divided by Column 4. 

The entries in Column 9 are the products of Column 6 and Column 3. 

The values presented in Column 6 have been calculated front the tem- 
perature spectra and the DO response to temperature, without use of the 
river discharge in the computation. It is highly significant that the mean 
of DO response at all the Potomac River Stations (1-8) is 0.15 ppm/ 
degree C., the exact rate of change in DO saturation concentration given 
for fresh water in standard tables (e.g., Standard Methods ). This rate 
of change applies at stations where the percent saturation of DO is low 
as well as where it is high. 

The values in Column 8 are the DO responses to river discharges, 
which include the advected DO as well as the change due to saturation 
concentration of DO, By subtracting the values in Column 9 from those 
in Column 8, a measure of the change in DO advected to the system is 
obtained. Station 0 shows a decrease in DO advected with increasing dis- 
charge, probably because the greater volume of water in the main stem 
reduces the amount of tidal exchange in this tributary. Since the tidal 
prism for the estuary above any point is a fixed volume, the amount of 
water entering and leaving an embayment such as the Anacostia River 
would be inversely proportional to the river discharge entering the sys- 
tem in a tidal cycle. 

A definite break in the magnitude of DO advected is noted between 
Stations 1 and 2 and the other stations on the main stem of the Potomac. 
The advection rates at Stations 1 and 2 indicate that the river water 
entering the estuary may be about 88 percent saturated with DO (based 
on Station 2, Column 9 divided by Column 10). An estimate of the tem- 
perature of the inflowing water may also be made. From the temperature 
response to river discharge at Stations 1 and 2 and the mean temperature 
values at these stations (28,0 and 28,1, respectively) the mean river tem- 
perature for the survey period must be about 26.6*C.; this Is certainly 
accurate within one degree. Temperature data from a water intake at 
Great Falls (about 30 river miles above the estuary) showed an average 
water temperature for this period of 26*C. The lower advected DO values 
for stations farther down in the estuary reflect the fad that some of the 
dissolved oxygen is used up before it can get to these stations, since advec- 
tive transport in an estuary is partly controlled by the reversible nature of 
tidal flow; a detailed analysis of these values, waste discharges, and tidal 
movement might, however, throw some light on the kinetics of waste assi- 
milation in this particular system. 

These results may be used to estimate the total pounds of oxygen 
entering the system each day in, and directly due to, the river flow. The 
purpose here has been merely to illustrate how the needed factors relating 
DO, temperature, and river discharge can be obtained by spectral analysis 
of the appropriate records. 



Example 5. Relationship cf River Flow and Tide to Stratification in an 

Estuary 

This example illustrates the use of spectral analysis in attacking the 
classic problem in estuarine hydrodynamics; that is, changes in stratifica- 
tion patterns with changes in river flow and tide range. 

Here also, a different approach to the use of spectral results is pre- 
sented. In the first example of this chapter a rationale was presented for 
choosing between the single estimate response and the overall response. 
In the remaining examples the single estimate response and the overall 
response have both been used, but with an emphasis on the single esti- 
mate response approach. That is, the progression of diumal effects was 
followed down the river or estuary, or the variations in semidiurnal tidal 
effects were considered at different points in an estuary. 

This present example deals entirely with overall response values. The 
reason for this is that the area of concern is the transfer of energy between 
the long-period-dominated river discharge and the short-period-dominated 
tidal effects. The magnitude of this energy transfer is what will force a 
change in stratification pattern; the shape of the spectrum itself will 
result from a variety of nonlinear interactions with which present theo- 
retical knowledge cannot cope. 

The overall response statistic compares one total energy spectrum with 
another and makes it possible to avoid the complications of nonlinear 
interactions in solving the present problem. 

The study area for this example is the harbor of Charleston, South 
Carolina, Figure 34; the data on which the analysis is based were obtained 
during six 5-day surveys at a variety of river discharges and tide ranges.* 
Each of the records analyzed consisted of 30 points at 4-hour sampling 
intervals; the analysis was run at 10 lags. 

Because of the shortness of the records, an analysis of variance was run 
on all combinations of records to be used in the spectral analysis. The 
significance of the F-ratios for each interaction was used as a guide for 
selecting environmental parameters and water quality parameters for spec- 
tral analysis. 

From the cross-spectral analysis, a coefficient for the overall response 
of chlorides to river discharge was computed. This coefficient indicates the 
magnitude of the effect of river discharge on chloride concentration and is 
expressed in units of gm/l/cfs. Table II contains the overall response 
coefficients calculated for each station and survey, together with the mean 
values of these coefficients for each survey. These coefficients include the 
effect of variations of all periods. 

Comparison of the results from surveys AA and AB showed that the 
means of surface and bottom responses of chloride concentrations t j river 
discharge were about equal for both surveys, 1.29 for AA and 1.35 for 
AB. On the surface the overall response was slightly higher for the spring 
tide condition, and on the bottom the response was slightly lower for the 



• The complete report from which this example h taken Is: ’Charleston Harbor 
Water Quality Study." Federal Water Pollution Control Administration, USD1, 
lune 1964. 
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spring tide condition. Mean surface-to-bottom response ratios of 0.555 
and 0.46, respectively, showed that the estuary was well stratified during 
both surveys. The mean for the two surveys was 0.51. These values com* 
pared well with the measured mean chloridt ratios for these surveys of 
0.57 (Table 12). 
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Table 11. Overall Response of Chloride Concentration to 
River Discharge, Charleston Harbor 









(mg/l/cf*) 






Survey 


AA 


AB 


B 


C 


D 


B 


Statioa 

IS 


1.292 


1.844 


.294 


.362 


.103 


.167 


IB 


1.935 


1.623 


.318 


.298 


.086 


.152 


2S 


1.513 


1.563 


.115 


.155 


.088 


.209 


2B 


1.475 


2.139 


.413 


.193 


.135 


.203 


3S 


.<427 


.500 


.088 


.157 


.052 


.084 


3B 


1.638 


1.852 


.411 


.328 


.111 


.133 


4S 


1.050 


.655 


.250 


.156 


.077 


.156 


4B 


2.186 


1.684 


.300 


.204 


.089 


.188 


5$ 


1.637 


.756 


.212 


.149 


.090 


.140 


SB 


2.317 


1.329 


.238 


.254 


.077 


.158 


6S 


.737 


.408 


.358 


.136 


.060 


.181 


6B 


1.509 


1.698 


.550 


.245 


.101 


.231 


7S 


1.008 


.758 


.124 


.089 


.094 


.166 


7B 


1.984 


1.762 


.292 


.214 


.137 


.171 


8S 


.633 


1.062 


.165 


.079 


.139 


.121 


8B 


2.451 


1.931 


.444 


.251 


.088 


.157 


9S 


.466 


,654 


.192 


.117 


.028 


.068 


9B 


.798 


2.073 


.358 


.172 


.053 


.088 


I OS 


.049 


.354 


— 


— 


— 


— 


10B 


.370 


2.434 


— 


— 


— 




ns 




— 


.030 


.117 


.049 


.096 


13B 


— 


— 


.456 


.384 


.110 


.206 


MeanS 


0.921 


0.855 


0,183 


0.152 


0.0812 


0.1436 


Mean B 


1.666 


1.853 


0.378 


0.254 


0.0974 


0.1646 


MeanS* B 


1.294 


1.354 


0.281 


0.203 


0.0893 


0.1541 


Ratios* B 


0.533 


0.461 


0.484 


0.598 


0,834 


0.872 



At high river discharges with the estuary strongly stratified the aver- 
age overall response of chlorides in the estuaiy was dependent upon the 
state of the tide. During spring tide conditions there was an increase in 
response in the surface layer which was related to more intense mixing 
at the interface created by higher shearing velocities. At the same time 
there would have been a stronger Intrusion of ocean water along the bot- 
tom, which meant less variation and consequently less response In the 
bottom layer. 

The response spectra for the other four surveys showed generally 
strong semidiurnal responses with the responses at each station becom- 
ing relatively stronger on the surface as the river discharge decreased. Con- 
currently, the overall response of chlorides to river discharge decreased 
for surveys B, C, and D, end then increased during survey E (made dur- 
ing spring tide conditions). These results reflected a decrease in the 
strength <5 the stratification with reduced river discharge and a consequent 
increased mixing between surface and bottom. The main ratios of chlo- 
rides as shown in Table 12 show the same type of progression. 

Surveys D and E show mean overall response ratios of 0.834 and 
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Table 12. Surface- to Bottom-Chloride Ratios, 
Charleston Harbor 



Surveys 

1965 


AA 

3/3-7 


AB 

3/23-27 


B 

6/21-25 


C 

7/19-23 


D 

8/16-20 


E 

9/20-24 


Slack 

Tide 

Runs 


Station 1 


.707 


.498 


.619 


.728 


.720 


.886 


.733 


Station 2 


.695 


.498 


.529 


.648 


.683 


.933 


.820 


Station 3 


.734 


.524 


.845 


.748 


.789 


.918 


.900 


Station 4 


.645 


.405 


.045 


.619 


.741 


.816 


— 


Station 5 


.710 


.343 


.520 


.721 


.653 


.868 


.927 


Station 6 


.558 


.212 


.552 


.462 


.621 


.808 


— 


Station 7 


Ml 


.333 


.497 


.580 


.715 


.794 


— 


Station 8 


.462 


.371 


.424 


.417 


.473 


.636 


— 


Station 9 


.838 


.509 


.678 


.748 


.830 


.885 


.915 


Station 10 


.890 


.764 


Discontinued 


• • • 


. • • 


. • . 


Station 13 


... 


... 


.371 


.420 


.473 


.645 


.6)2 


M;an Ratio 
of 1 thru 9 


.670 


.410 


.572 


.630 


.692 


.838 


.784 


Mean Ratios 


AH Stations 


.692 


.445 


.552 


.609 


.670 


.819 




Average River 
Discharge cfs 


27,600 


27,700 


26,300 


19,300 


16,200 


13,500 


— 



0.872, respectively, indicating that the surface and bottom layers reacted 
nearly the same to river discharge for both surveys. 

The striking differences in overall response between surveys D and B 
required some examination. Table 13 shows the mean tide height range 
for each survey together with an analysis of the mean overall responses. 
The tidal height range difference between surveys D and E was about the 
same as that between surveys AA and AB. The surface and bottom re- 
sponses for surveys AA and AB were quite different from those during 
surveys D and h. However, the difference in surface response between 
both pairs of surveys was practically identical, showing that the .surface 
layer was being affected similarly by some constant difference between 
surveys comprising each pair. The only similarity in environmental fac- 
tors was the change in tidal range. 

One additional result is important; between surveys D and E the change 
in bottom response was the same as the change in surface response, 0.062 
on the surface and 0.067 on the bottom. This comparison shows that 
the surface and bottom overall responses for surveys D and E were very 
nearly the same and that the changes of overall response with change 
in tide height were the same in both surface and bottom records. The net 
result of this comparison of chloride responses was to demonstrate that 
surface and bottom layers in surveys D and E behaved so nearly alike 
that, for all practical purposes, the chloride distribution is acting as if 
estuary stratification was virtually nonexistent. 

The cross-spectral analysis further demonstrated that, below discharges 



Table 13. Analysis of Chloride Responses to Tide Ranges, 
Charleston Harbor 













Differences 






Mean Tide 


Surface 


Bottom 


Tide 


Surface 


Bottom 




Range 


Response 


Response 


Range 


Response 


Response 


Survey 


(ft.) 


mg/ 1 /cfs 


mg/l/cfs 


(ft) 


mg/l/cfs 


mg/ 1 /cfs 


AA 


5.22 


0.921 


1.666 


1.14 


.066 


0.187 


AB 


4.08 


0.855 


1.853 








B 


4.09 


0.183 


0.378 








O 


4,11 


0.152 


0.254 








D 


4.42 


0.0812 


0.0974 


1.34 


.062 


0.067 


E 


5.76 


0.1436 


0,1646 









of about 16,000 cfs, the Charleston Harbor system, considered as a sin- 
gle unit, behaved in its chloride distribution as if it were vertically un- 
stratified. 

For analyzing the chl< .ide — river discharge relationship at flow rates 
below the range observed during the intensive survey program, the chlo- 
ride response factors used were the average of those found during surveys 
D and E (0.834 and 0.872, respectively) because conditions dose to 
vertical mixing were observed during these surveys. This would corre- 
spond to a tide range of 5.0 feet, which Is near the mean tide range of 5.2 
feet for the primary reference tide in Charleston Harbor. 
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